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Abstract: We derive a number of local identities of arbitrary rank involving Jacobi elliptic functions and 
use them to obtain several new results. First, we present an alternative, simpler derivation of the cyclic 
identities discovered by us recently, along with an extension to several new cyclic identities of arbitrary 
rank. Second, we obtain a generalization to cyclic identities in which successive terms have a multiplicative 
phase factor exp(2i7r/s), where s is any integer. Third, we systematize the local identities by deriving four 
local "master identities" analogous to the master identities for the cyclic sums discussed by us previously. 
Fourth, we point out that many of the local identities can be thought of as exact discretizations of standard 
nonlinear differential equations satisfied by the Jacobian elliptic functions. Finally, we obtain explicit 
answers for a number of definite integrals and simpler forms for several indefinite integrals involving Jacobi 
elliptic functions. 
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1 Introduction 



In a recent paper [I], (henceforth referred to as I), we have given many new mathematical identities involving 
the Jacobi elliptic functions sn(x,m), cn(x,m), dn(x,m), where m is the elliptic modulus parameter 
(0 < m < 1). The functions sn(x, m), cn(x,m), dn(x,m) are doubly periodic functions with periods 
(4K(m),i2K'(m)), (4K(m),iAK'(m)), (2K (m),i4if' (m)), respectively HE]. Here, K(m) = J* /2 M[l- 
m sin 2 denotes the complete elliptic integral of the first kind, and K'(m) = K(l — m). The m = 

limit gives K(0) = it/2 and trigonometric functions: sn(x,0) = sinx, cn(x,0) = cosx, dn(x,0) = 1. The 
m — > 1 limit gives K(l) — ► oo and hyperbolic functions: sn(x, 1) — ► tanhx, cn(x, 1) — ► sechx, dn(x, 1) — ► 
sechx. For simplicity, from now on we will not explicitly display the modulus parameter m as an argument 
of the Jacobi elliptic functions. 

The identities discussed in ref. I are all cyclic with the arguments of the Jacobi functions in successive 
terms separated by either 2K(m)/p or 4K(m)/p, where p is an integer. Each p-point identity of rank R 
involves on its left hand side a cyclic homogeneous polynomial in Jacobi elliptic functions of degree R with 
p equally spaced arguments. The separation is 2K(m)/p or 4K~(m)/p depending on whether the period 
of any term on the left hand side is 2K(m) or 4K(m). In another recent publication [1] (referred to as 
II), we presented rigorous mathematical proofs valid for arbitrary p and R even though, for simplicity, we 
only presented identities of low rank. In ref. II, we classified the identities into four types, each with its 
own "master identity" which we proved using a combination of the Poisson summation formula and the 
special properties of elliptic functions [IJ[3]. We also provided a rigorous derivation of cyclic identities with 
successive terms having alternating signs. 

In this paper, we provide several generalizations of the results discussed in refs. I and II. Here, our 
approach is different and involves the use of "local" identities which focus on just any one term in a 
cyclic identity. This term involves a product of Jacobi elliptic functions and is expressed via the local 
identity as the sum of many terms of lower rank. The purpose of this paper is to derive and make use of a 
number of local identities for Jacobi elliptic functions. These local identities form the building blocks for 
cyclic as well as much more general identities. For instance, adding p local identities with equally spaced 
arguments permits us to re-derive cyclic identities. More generally, taking p local identities with a phase 
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(-1)0-1) = e (j-l)in and 

summing over the index j gives previously derived identities in which successive 
terms have alternative signs. Finally, as discussed below, the generalization to taking p local identities with 
an even more general phase exp(2i(j — 1)tt/s), where s is any integer and summing over the index j yields 
interesting new identities in which successive terms have different weights. Also, while in principle we were 
able to prove the general form for identities of arbitrary rank in refs. I and II, in practice it was very 
difficult to obtain the explicit coefficients in these identities. The use of local identities permits evaluation 
of these coefficients. As a byproduct, a number of definite integrals involving Jacobi elliptic functions can 
be explicitly evaluated and a number of indefinite integrals can be expressed in simpler form. Finally, 
we show that some special linear combinations of various identities mentioned above have a particularly 
simple right hand side, and we give several illustrative examples. 

To clarify the above ideas about the approach used in this paper, consider as an example, one specific 
basic local identity derived here: 

dn 2 (y)dn(y + a) = — cs 2 (a)dn(y + a) + ds(a)ns(a)dn(y) — mcs(a)cn(y)sn(y) . (1) 

Choosing y = x + (j — l)2K(m)/p with j = 1,2, ...,p actually corresponds to p identities, one for each 
value of the integer j. Taking a = r2K(m)/p, where r is an integer which is less than p and coprime to it, 
and summing over j yields the cyclic identity 



i=i i=i 



-dj - mc8(a)8jCj 



(2) 

where the coefficient A is given by A = 2[ds(a)ns(a) — cs 2 (a)] and we have used the notation 

dj = dn(x + {j-l)2K(m)/p,m), sj = sn(x + (j - l)2K(m )/p, m), cj = cn(x+ (j - l)2K(m)/p, m) . (3) 

Similar manipulations using a = —r2K(m ) jp yield cyclic identities for expressions like Y%=i dj [rfj+r±dj-r] • 
The result with the negative sign is new and will be discussed later in this paper. The result with the 
positive sign is 

v p 

d] [dj+r + dj- r ] =Aj2d 3 , (4) 

3=1 3=1 

and is one of the cyclic identities derived in ref. II by more complicated techniques. 
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Further, while it was clear from refs. I and II that the identities of arbitrary (odd) rank which are 
generalizations of eq. (£Q) must have the structure 

J2 df [d j+r + dj-r] = A l Y J df~ 1 +... + A n J2 dj , (5) 

j=i j=i j=i 

we were unable to obtain the coefficients Ai,...,A n . Here, we will obtain explicit expressions for the 
coefficients. In ref. II, we were able to obtain the analogue of identity (j3J) with alternating signs given by 



p p OK 
5^(-iy _1 t§ [d j+r + dj-r] = A^-l)' -1 ^ , A = 2 ds(a)ns(a) + cs 2 (a) , a = . (6) 

j=l j=i P 

Here, we will obtain identities with more general weights like 

P ._ 2iir 

^2uj j 1 d"j[dj +r + dj-r] , uj = exp( ), (7) 

3=1 S 

where to is the s th root of unity, with s being any integer (< p) and p being mod s. Finally, in ref. II we 

had obtained MI-II (class II master identity, also see Sec. 4 below) identities like 

P P n ( f 2K \ r2 K 

YdM = -2cs 2 {a)Y d 2 , + ^- / dn 2 (t)dn 2 (t + a)dt + 4Ecs 2 {a) , a= , (8) 

JZi JZi J yJo J p 

where E is the complete elliptic integral of second kind : Z\- The approach in this paper will permit an 
evaluation of the definite integral on the right hand side. 

The plan of this paper is as follows. In Sec. 2, we state several local identities and indicate how they 
are derived. It may be noted here that each identity has an integer label R indicating the rank of the 
identity, i.e. the left hand side of the identity is a homogeneous polynomial of degree R. We also show here 
that linear combinations of cyclic identities often yield simpler results. In Sec. 3, we use the local identities 
of rank 2, 3, 4 recursively to obtain local identities of arbitrary odd and even rank, using which one can 
immediately obtain the corresponding cyclic identities with weight factors uj. In Sec. 4, we provide a unified 
framework for the local identities by deriving four master local identities, from which all the identities can 
be derived in an alternative manner without using addition formulas. In Sec. 5, we concentrate on those 
identities of ref. II in which one of the terms on the right hand side is a definite integral (which we were 
previously unable to evaluate). Using our local identities, we show that one can obtain cyclic identities 
where all the terms on the right hand side are now explicitly known. In fact, we show that by starting 



from any given local identity, the indefinite integral of the left hand side of this identity can be analytically 
obtained in terms of well known integrals of Jacobi elliptic functions and indefinite elliptic integrals of 
the first, second and third kind j3 El- We would like to re-emphasize that most of these integrals do 
not seem to be known in the literature. In Sec. 6, we discuss continuum limits of the local and cyclic 
identities, showing that these degenerate to standard differential equations or integral formulas. Sec. 7 
contains conclusions and a discussion of some open problems. All local identities of ranks 2, 3 and 4 are 
presented in Appendices A, B and C respectively. A few local identities of rank 5 and arbitrary rank are 
presented in Appendices D and E respectively. Several simple results obtained by taking suitable linear 
combinations of cyclic identities are given in Appendix F. Many new definite and indefinite integrals are 
given in Appendices G and H respectively. 



2 The Basic Local Identities 



In this section we shall obtain several basic local identities. These identities are easily derived using the 
well-known addition formulas for the sn, cn, dn functions [HE]: 

dn(a)dn(6) — mcn(a)cn(6)sn(a)sn(6) 



dn(a + b) = 
cn(a + b) 
sn(a + b) 



1 — msn 2 (a)sn 2 (6) 

cn(a)cn(6) — dn(a)sn(a)dn(6)sn(6) 

1 — msn 2 (a)sn 2 (5) 
sn(a)cn(6)dn(6) + cn(a)dn(a)sn(6) 
1 — msn 2 (a)sn 2 (6) 

We shall also use the addition formula for the Jacobi zeta function given by 

Z(a + b) = Z(a) + Z(6) — msn(a)sn(6)sn(a + b) . 
One of the simplest, local, rank two identities is 

dn(x)dn(a; + a) = dn(a) + cs(a)[Z(x + a) — Z(x) — Z(a) 



(9) 
(10) 
(11) 

(12) 
(13) 



which is easily proved by algebraic simplification after using the addition formulas © to (|12|) . The power 
of this local identity can be appreciated by the fact that we can immediately derive the cyclic identity 



J2djd j+1 = p 



,2K X ,2K, ,2K S 
dn _ -cs — )Z — 
p p p 



(14) 



which was obtained in ref. II. This is done by writing local identities like (j!3|) with x being replaced by 
x + a, x + 2a, . . . , x + {p — l)a and choosing a = 2K/p. On adding these identities and noting that 
dn(x) has a period 2K, we then immediately obtain the cyclic identity (|14|) , Here p denotes the number of 
subdivisions of the period at which Jacobi elliptic functions dn(x) are evaluated. A generalization of this 
identity to r th neighbours is immediate, i.e. on choosing a = r2K/p (where r is coprime to and less than 
p), we obtain the more general identity 

p r2K 

2_, djdj+r = p[dn(a) — cs(a)Z(a)] , a = , (15) 

j= i ' P 

which was also obtained in ref. II. 

We can immediately obtain a local identity for dn(x)dn(x — a) by changing a to —a and recognizing 

the fact that while cn(a),dn(a) are even functions of a, the functions sn(a),Z(a) are odd: 

dn(x)dn(x — a) = dn(o) — cs(a)[Z(x — a) — Z(x) + Z(a)] . (16) 

Adding and subtracting eqs. (fT3|) and (fTO)) yields alternative simple expressions: 

dn(x)[dn(j; + a) + dn(x - a)] = 2dn(a) + cs(a)[Z(x + a) - Z(x - a) - 2Z(a)] , (17) 

and 

dn(x)[dn(x + a) - dn(x - a)} = cs(a)[Z(sc + a) + Z(x - a) - 2Z(x)} . (18) 

If we now consider the local identities analogous to (fT3)l with x being replaced by x + a, x + 2a,..., x + 
(p — l)a, multiply them in turn by u, uj 2 ,..., u> p ~ 2 , ujP" 1 respectively and add to the local identity (|13|l. 
then we obtain the remarkable identity 

y^w J l djdj +r = p[dn(a) — cs(a)Z(a)]5 s i — I 1 I cs(a) 1 Z J - , (19) 

J l \ j=1 

where a = r2K/p. The phase a; is as given by eq. with s < p and p being mod s. For the special case 

s = 1 we recover the cyclic identity ()15|) with all terms on the left hand side having positive signs. For 

s = 2, eq. (|19|) gives the cyclic identity with terms having alternating signs as obtained in ref. II. Thus 

the local identities are very basic in the sense that once they are known, then the corresponding cyclic 
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identities with and without arbitrary weight uj are immediately obtainable. It is worth emphasizing here 
that the cyclic identities with arbitrary weight are new. 

Proceeding in the same way, we have derived all possible local identities of rank two, three and four. 
They are given in Appendices A, B and C respectively. Some examples of local identities of rank 5 are 
given in Appendix D. By following the procedure explained above, in each case it is easy to obtain the 
corresponding cyclic identities with weights co. 

One advantage of the local identities approach is that for the MI-II type of cyclic identities, the right 

hand side is explicitly known. In this context, it is worth mentioning that in ref. II (also see [Hj) we 

had obtained several MI-II cyclic identities in which one of the terms on the right hand side is a definite 

integral. For example, one of the cyclic MI-II identities obtained in ref. II is given by eq. (jSJ). However, 

if we take the local identity (|166|) given in Appendix C, and use the procedure described above, we find a 

simpler, more elegant form for this cyclic identity 

p p r2 R 

d 2 d 2 j+r = -2cs 2 (a) ^ d 2 + p[cs 2 (a) + ds 2 (a) - 2cs(a)ds(a)ns(a)Z(a)] , a = . (20) 

j= i j= i P 

Proceeding in this way and using various local identities obtained in this paper, the corresponding MI-II 
cyclic identities are easily written where the constant on the right hand side is now explicitly known and 
is not just formally expressed as an unevaluated definite integral. 
Let us now consider the following cyclic identity of rank two 

mcn(x)[sn(x + a) — sn(x — a)] = 2ns(a)dn(x) — ds(a)[dn(x + a) + dn(x — a)] , (21) 

which is easily derived by using eqs. © to (|12[), ^From here we easily obtain the following cyclic identity 
with weighted terms 

P r 2tt 1 p 

^maj J l Cj[sj +r — Sj- r ] = 2 ns(a) — cos( — )ds(a) 1 dj , (22) 



where a = r2K/p. Some examples of new cyclic identities (with weighted terms) of rank 3 and 4 are: 

p 

m oji~ 1 dj [cj +r Sj +r — Cj_ r Sj_ r ] = 2p[cs(a) — ds(a)ns(a)Z(a)]5 s i 
3=1 

2ir v 2vr p 

-2isin(— )ds(a)ns(a) -2cos(— )cs(a) ^V^df , (23) 

S 3=1 S 3=1 



^u;- 7 l dj[cj +T dj +r — Cj- r dj- r } = 2cs(6) cos( — ) Sjdj — 2« sin( — )ds(6)ns(6) cj- 7 Cj , (24) 



3=1 

,2tt p 
s 



27r > 

-2i sin( — )cs(6)ns(6) ^ uj^ l Sjdj 
s . 

3=1 



-2m cos(— )ds(6) + 2ds ( b ) 

3=1 



2vr, 



(m + cs (6)) cos( — ) — cs(6)ns(6) 



p rAK 



It may be noted that the identities and (|23|) are of type MI-IV and hence Cj = cn(x + (j — l)4K/p, m). 
Proceeding in the same way, corresponding to most of the cyclic identities discussed in I and II, we can 
obtain the corresponding cyclic identities with weighted terms. For illustration purposes, a few such cyclic 
identities are presented in Appendix F. 



3 Local Identities of Arbitrary Rank 

Now that we have obtained local and hence cyclic identities of low rank, the obvious question to ask is 
whether one can generalize and obtain corresponding local (and hence also cyclic) identities of arbitrary 
rank. In this section, we show that this is indeed possible. In particular, we show that corresponding to 
each local low rank identity, we can obtain a local identity of arbitrary rank in which all the coefficients 
are explicitly known. 

As an illustration, let us start from the local identity 

dn 2 (x)[dn(x + a) + dn(x - a)} = Adn(x) + B[dn(x + a) + dn(x - a)] , (26) 

where A, B are constants 

A = 2ds(a)ns(a) , B = -cs 2 (a) . (27) 

This identity can be easily derived using the addition formula ©. On repeatedly multiplying both sides 
of identity (|26|) by dn 2 (x) and using eq. (|26[) we obtain the following local identity of arbitrary odd rank 

n 

dn 2n (x)[dn(x + a) + dn(x -a)] = A^ B k - l dn 2( - n -^ +1 (x) + B n [dn{x + a) + dn(x - a)] , (28) 

k=l 
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where the constants A, B are as given by eq. (|27|). Using the procedure described in Sec. 2, the corre- 
sponding cyclic identities with and without arbitrary weight uj are immediately obtained: 

P P n P 

Y / dr[d j+ r + d j - r ]=Aj2 J2 Bk - ld f n ~ k)+1 + 2Bn Y, d l > ( 29 ) 
j=l j=l k=l j=l 

P p n 9 V 

J2 ^-^[dj+r + dj-r] =^EE u^B*" 1 ^ n ~* )+1 + 2B" cos(— ) ^' l dj , (30) 

j=l j=l k=l S j=l 

where A, B and uj are as given by eqs. (|27jl and (JZJ) respectively while a = 2rK/p and p = mod s. Note 
that the identity without any weight (i.e. (|29|l) is obtained from ()3fl)) in the limit s = 1 while for s = 2 
we obtain the identity with alternate sign. It is worth emphasizing that in the identities (|28jl to (|30|) of 
arbitrary odd rank, all the coefficients are explicitly known. 

In order to obtain the corresponding local identity of arbitrary even rank, we start from the identity 
(|28|) and multiply both sides of it by dn(x) and use the local identity (fT7|) to obtain 

n 

dn 2n+1 (a;)[dn( a ;+a)+dn(2;-a)] = A^ B k - l dn 2 ^- k+1 \x)+2B n dn(a)+B n cs(a)[Z(x+a)-Z(x-a)-2Z(a)]. 

k=l 

(31) 

The corresponding cyclic identities with and without arbitrary weight are then immediately obtained by 
following the steps as above. 

Proceeding in the same way, but starting from the identity 

dn 2 (x)[dn(x + a) — dn(x — a)] = Dcn(x)sn(x) + B[dn(x + a) — dn(x — a)] , D = — 2mcs(a) , (32) 

multiplying recursively by dn 2 (x) and using the identity (|3*2*|) , we obtain the following identities of arbitrary 
odd and even rank 

n 

dn 2n (x)[dn{x + a) - dn(x - a)] = B k - 1 cn(x)sn(x)dn^ n ~ k Hx) + B n [dn(x + a) - dn(x - a)] , (33) 

k=l 

n 

dn 2n+1 (x)[dn(x+a)-dn(x-a)] = B^ 1 cn(x)sn(x) 

dn 2(n-fc)+i ^ +jB n cs ( a ) [z( x +a) +Z(x-a) - 2Z(x)] , 

k=l 

(34) 

where B = — cs 2 (a) and D = — 2mcs(a). The corresponding cyclic identities are then immediately written 
down. Further, by adding the two identities ()28|) and H33j) . we obtain the basic local cyclic identity of any 
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odd rank: 

D n A n 

dn 2n (x)dn(x + a) = - £ B k - 1 cn(x)sn(x)dn 2 ^- k \x) + £ n dn(x + a) + - ^ B^dn 2 ^^^) , (35) 

2 fe=l 2 fe=l 

where A,B,D are given by eqs. I|27[l and (|32|) . It is worth pointing out that using this identity, we can 

immediately write down the local identity for the combination dn(x)dn 2n (x + a). This is done by replacing 

x by x — a followed by changing a to —a in eq. (|35|). In this way we obtain 

A n D " 
dn(x)dn 2n (x+a) = -^B^W^^+^x+a) ^ fc_lcn ( x + a ) sn ( x + a ) dn2(n_fe) (a;+a)+ J B n dn(x). (36) 

2 k=i 2 k=i 
Now dn(x)dn 2n (x — a) can be immediately obtained from here by replacing a by —a. We find that 

v 



£ df [d j+r ± d^-r] =±£ d i ± d f-r] . ( 37 ) 



^(-i)i-^f [d i+r ± = T E(-ir^ ± 4-] • (38) 

j=l 3=1 

In the next section we shall see that similar relations are in fact true in general for any such combinations 
of Jacobi elliptic functions. 

Proceeding in the same way, by starting from each of the lower rank identities given in Appendices A, 
B, C we can write down the corresponding identities of arbitrary even as well as odd rank. Some illustrative 
examples with arbitrary even as well as odd powers of dn(x) (or sn(x) or cn(x)) are given in Appendix E. 



4 Master Local Identities 

In ref. II, we derived four master identities from which all the cyclic identities could be derived as special 
cases. In this section we show how a similar procedure works at the level of the local identities, thereby 
systematizing the identities and providing a unified framework for them. Besides, rather than using the 
addition formulas for Jacobi elliptic functions, the master identities (MI) provide an alternative way to 
derive the local, and hence also cyclic, identities. We may note here the differences that arise in the two 
approaches: addition formulas do not lead to unevaluated constants in the form of integrals on the right 
hand side, while the MI, in particular one of the four classes of MI do; on the other hand MI reduces the 
right hand side maximally to standard forms, while the addition formulas approach needs considerable 
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algebraic manipulation to attain the simplest final form. In any case, the two approaches are of course 
compatible and either of them may be used. In this section and in Sec. 6, we will use z instead of x as the 
variable to emphasize that there is no restriction to the real numbers. 

The classification in ref. I and in the above sections has been in terms of the polynomial order of the 
elliptic functions appearing in the left hand side, called the rank of the identities. In contrast, the master 
identities first use symmetries to identify four classes. Within each class, the identities are characterized 
by a number which is the highest order of the singularities in the fundamental domain of the left hand side. 
Thus, the analytic structure of the functions appearing on the left hand side of any identity determines 
the constants and the form of the functions that appear on the right hand side. It is then quite easy to 
write symbolic manipulation programs that turn any given form of the left hand side into an appropriate 
local identity. 

We first recall essential details of the analytic properties of the Jacobi elliptic functions _6J . The function 
dn(z) is an even elliptic function of order two; there are two simple poles inside the period parallelogram 
(0, 2K, 2K + 4iK' ,4iK') situated at iK' and 3iK' with residues of — i and i respectively. The function sn(z) 
is an odd elliptic function of order two; with two simple poles situated at iK' and iK' + 2K, with residues 
l/^/m and — l/^/m", inside the fundamental period parallelogram (0, AK, 4K + 2iK', 2iK'). The function 
cn(z) is an even elliptic function of order two; with two simple poles situated at iK' and 2K + iK' , with 
residues —i/y/rK and i/^/m, inside the fundamental parallelogram (— 2K, 2K, AK + 2iK', 2iK'). We note 
that the lattice of the poles in the complex plane is identical for all these three functions. However these 
functions have the following important distinguishing properties that we will use below: dn(z + 2iK') = 
-dn(z), sn(z + 2iK') = sn(», cn(z + 2iK') = -cn(z), dn(z + 2K) = dn(z), sn(z + 2K) = -sn(z), 
cn(z + 2K) = -cn(z). 

The symmetry and periodicity properties put together allow us to concentrate on the region (0, 2K, 2K+ 
2iK' , 2iK') uniformly for all the functions, and consider only one simple pole at iK' . We supplement these 
possible symmetries with one additional one, for which we describe the properties of the elliptic function 
dn 2 (z). Equivalently one may choose cn 2 (z) or sn 2 (z). The function dn 2 (z) has the fundamental domain 
(0, 2K, 2K + 2iK' ,2iK') and consequently is completely periodic with respect to translations of 2K and 
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2iK' . It is also of order two, with one double pole at iK' , with a residue of 0. Thus we classify functions 
f(z) constructed from the Jacobian elliptic functions into four symmetry classes. We define the quantities 
P,Q by: 

f(z + 2iK') = (-l) p f(z), f(z + 2K) = (-lff(z), P,Q = 0,1. (39) 

We denote the four possibilities by (—,+), (+,+), {+,— ) and (—,—), where the first sign refers to 
the sign of (— l) p and the second to that of (—1)^. We note that as far as periodicity is concerned these 
functions are identical to dn(z), dn 2 (z), sn(z) and cn(z) respectively. We also note that repeated differ- 
entiation does not change the symmetry class to which the functions belong, while this creates functions 
with arbitrarily high order of poles. This then allows us to tailor suitable combinations of derivatives of 
these four functions such that not only the periodicity, but also the singular parts match with the given 
function f(z). Thus the difference between the function f(z) and the tailored combination is an elliptic 
function with no poles anywhere including at infinity. We then use Liouville's theorem that states that 
if an analytic function has no pole anywhere including at infinity then it must be a constant, and then 
explicitly show that the constant is zero in all cases except the second, where it can be evaluated as a 
definite integral. 

4.1 Master Identities of Types I, III and IV 

Let f(z) be an elliptic function with the symmetry properties corresponding to P = 1, Q = 0, (type I) and 
having n p poles at positions a r (r = 1, . . . , n p ) within the region (0, 2K, 2K + 2iK', 2iK') which we will 
call ABCD. Let the principal part around the pole a r be 

Y. i — (40) 

We note that the principal part of dn(z) around the pole iK' is 



—i 



{z - iK') 

Therefore if we consider the function g(z): 



(41) 



U P L r ., n L-l , , ,L.-1 



i(-l) r ( r )d lr 1 dn(z). /A . 

r=l l r =l v ' 
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this has identical poles as f(z) and at these poles also has identical principal parts. Due to the symmetry 
requirements, the functions f(z) and g(z) also have identical periods and hence they by Liouville's theorem 
they can differ utmost by a constant that is independent of z. However integrating both these functions 
from to 4iK' , we see from the antisymmetry that these must vanish, implying that the constant is zero; 
and hence f(z) = g(z). This is our "master" local identity of type I; often the evaluated function g(z) is 
of a simpler form than f{z). 

Consider as an illustration the identity that results when f{z) = dn 2 (z)[dn(z + a) + dn(z — a)]. This 
function has three poles, within ABCD at a\ = iK' , a<i = —a + iK' , and = a + iK' . This function has 
P = 1, Q = 0, and hence is of type I. At a\ = iK' the principal part is 

-2ids(a)ns(a)/(z - iK'). 

We note that although dn 2 (z) has a double pole at iK' it gets "softened" by one, because dn(z+a)+dn(z— a) 
has a zero at iK' for all a. This is the reason why we expect that the RHS of such identities are simpler 
than the LHS. Thus = — 2i ds(a)ns(a) and L\ = 1. Similarly we get: = ics 2 (a), L2 = 1 and 
a>i = ics 2 (a), L3 = 1. Hence this master identity yields the already stated result in Eq. Q26j) . which was 
alternatively derived using addition formulas. 

We note that in the case of cyclic identities further simplification occurs and the a at the various poles 
can be summed up, while in the case of local identities they are left as they are. We note that the structure 
of the LHS yielded to simplification because roughly a zero cancelled a pole. We can look at the "parts" 
of this identity where this does not happen fully. Thus when we take f(z) = dn 2 (z)dn(.z — a) we get the 
identity: 

dn 2 (z)dn(z — a) = ds(a)ns(a) dn(z) + m cs(a) cn(z)sn(z) — cs 2 (a)dn(2: — a). (43) 

We note that the rank of the RHS is one less than that of the LHS, and further reduction by one occurs 
when a is changed to —a and the two identities for dn 2 (z)dn(z + a) and dn 2 (z)dn(z — a) are added which 
results in the already quoted identity of Eq. (|26|) . 

Similarly we derive the master identity for functions belonging to type III (P = 0, Q = 1): 

f( Z ) = EE V-l), -V^^U-ar^, (44) 
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and type IV (P = 1,Q = 1): 

* £ i^(-l)^ 1 (r) jr-l cn (,) 

/(*) = E E (fr _ 1}1 < -sczrU w (45) 

r=l Z r =l v ' ' 

That the constant is zero in the case of type III and type IV identities can be seen by integrating both 
sides from to AK. It may be noted that symbolic manipulation packages that calculate series expansions 
can be effectively used to generate these identities. 

4.2 Master Identity of Type II 

This last type of identity deserves special mention; firstly the function archetype is dn 2 (z) which has a 
double pole at iK', secondly it leads to identities with non-zero constants, and lastly the Jacobian zeta 
function appears in an essential way. A function belonging to this type is periodic with periods 2K and 
2iK' (hence P = 0, Q = 0). Its principal part around iK' is 

-l/(z-iK') 2 . 

Thus we write the master identity in this CclSG clS 

/(*) = C + LL n _l)| a t dz l r -2 \z~a r +iK>, (46) 
r=l l r =l r >' 

where C is a constant. Note that the derivative order starts from —1, which should be interpreted as 
an integral. Functions of this type can have also simple poles and therefore the function dn 2 (z) and its 
derivatives are not sufficient to construct these. If we include its integral the master identity is complete, 
therefore we note the standard result [S] which can be taken to be the definition of the Jacobian zeta 
function Z{z): 

Z(z) 



o 



dn 2 (n) - | 



du = E(z) - ^z, (47) 



where E{z) is the incomplete elliptic integral of the second kind and E and K are the complete elliptic 
integrals of the second and first kinds respectively. Thus Z(z) is closely related to the incomplete elliptic 
integral of the second kind, and is periodic with a period of 2K, but is not elliptic. It is however almost 
elliptic due to the identity |SJ: Z(z + 2iK') = Z(z) — iir/K. If we write the l r = 1 part of this master 
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equation, which has the only part with the Jacobian zeta function, it is 

n p 

X^ai r) Z(z-a r + iK'). (48) 

r=l 

We note that this is an elliptic function with the correct periods of 2K and 2iK', due to the fact that: 

rip 

E«! r) = °- ( 49 ) 

r=l 

This is the sum of the residues of the function at all the poles in ABCD. Making use of the double 
periodicity of f(z) we find the integral around ABCD vanishes, hence from Cauchy's theorem it follows 
that the sum of the residues must also vanish, hence proving the above. Thus we are justified in using the 
zeta function in this type of master identity even if it is not elliptic: it will always appear in combinations 
that are elliptic functions. 

Integrating both sides of the master identity from to 2K we get to evaluate the constant C: 

r 2K 



1 f 2K 7o£ 



C 2K 

where 72 = Y%=i a 2^- We can therefore evaluate the identity at some convenient z where there is no 
singularity (for instance perhaps z = 0) and then make use of the above to evaluate definite integrals. 
Assuming for instance that there is no pole at z = we may write: 

1 f 2 %/u f(M E (-l)'- 1 (,) ^- 2 dn 2 (z) , 

— y o mdz = /(o) - 72 - - eEttttt)!^ ^-2 I*'-* ( 51 ) 

We note that this integral cannot be evaluated by a direct application of Cauchy's theorem due to the 
vanishing of both the contour integral around ABCD and its residue. This is an useful way of integrating 
many functions. 

We also point out here something that is of relevance to the cyclic identities: 

f>(^)[^+i) =±t%)bM±sMl, (52) 

3=1 3=1 

E(-l)Mz,)[M^+i) ± ^-1)] = T^t{-iyK Zj )[g{z 1+1 ) ±5(^-1)], (53) 
j=i 3=1 

where h(z) and g(z) are combinations of Jacobian elliptic functions as above. These relate ordinary and 
alternating sums under an interchange of h and g and follow from quite general considerations related 
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only to the periodicity of h and g. The RHS amounts to a rewriting of the LHS if we recall that either 
h(zp+i) = h(z%) and g(z p+ i) = g{z\) or h(z p +i) = —h(z%) and g{z v+ {) = —g(z\) as the whole function 
g(z)[h(z + T/p) + h(z — T/p)] is periodic with period T . 

5 Evaluation of Several Elliptic Integrals 

In ref. II, we obtained several cyclic identities of type MI-II in which the right hand side contained a definite 
integral involving products of Jacobi elliptic functions. These integrals are not available in standard tables 
of integrals [3J|5]. We now show that using local identities we can explicitly evaluate many such definite 
integrals. 

As an illustration, we start from local identity (|166j) . On integrating both sides with respect to x over 
an interval [0, 2K] yields the definite integral 

/ dn 2 (x)dn 2 + a) dx = -4Ecs 2 (a) + 2K[cs 2 (a) + ds 2 (a) - 2cs(a)ds(a)ns(a)Z(a)] . (54) 
Jo 

It may be noted that a is here any non-zero constant. Using this value of the integral in the MI-II cyclic 
identity (jHJ) that we obtained in ref. II and choosing a = immediately yields the cyclic identity (|20(1 
which we had directly obtained from the local identity. 

The other definite integrals which we are now able to evaluate are related to cyclic identities containing 
an even number of dn or sn or cn. For example, in ref. II, we simply stated the identity 

I P l r2K 

- djdj +r dj +s dj +t = A = —— / dn(x)dn(x + a)dn(x + a')dn(x + a") dx , (55) 
P ~[ 2K Jo 

but were unable to evaluate the integral and hence find A. Here, a = 2 ^,a' = 2 ^-,a" = 2 i^-. We now 
show that using the local identities, we can in fact evaluate this definite integral even for arbitrary but 
unequal a, a', a". To this purpose, we start from the local identity (|115|) . Integrating both sides over the 
interval [0, 2K] yields 

1 f 2K 

—— / dn(x)dn(x + a)dn(x + a)dn(x + a") dx = dn(a)dn(a')dn(a") + cs(a)cs(a' — a)cs(a" — a)Z(a) 
2K Jo 

- cs(a')cs(a' - a)cs(a" - a')Z(a') + cs(a")cs(a" - a)cs(a" - a')Z(a"), (56) 
which is eq. ()3U1(1 . Note that the special case a = , a' = , a" = yields the cyclic identity (|53|) . 
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Finally, there are some MI-II cyclic identities and hence definite integrals which were not even discussed 
in ref. II. For example, consider the local identity (|1U2|) from which by following the method explained in 
Sec. 2 we can deduce the cyclic identity 

1 P 1 f 2K 

- mdjSj +r Cj +s = — — / mdn(x)sn(a; + a)cn(x + a') dx, (57) 
Pj\ 2A Jo 

where a = ^y- , a 1 = . By using the local identity (|1U2() , we can in fact obtain this integral for arbitrary 
but unequal values of a, a'. In particular, on integrating both sides of eq. (|1U2|) over the interval [0,2iT] 
yields 

l r 2K 

— — / mdn(a;)sn(a;+a)cn(a;+a ) dx = — ds(a— a )[dn(o)— cs(a)Z(a)]+ns(a— a )[dn(o )— cs(a )Z(a')]. (58) 
2iv Jo 

In the special case when a = 2 ^ , a' = , one recovers the cyclic identity (|57]> . 

Proceeding in the same way, we have been able to obtain expressions for all the definite integrals which 
appear in the cyclic identities of the type MI-II (and which to the best of our knowledge are not known 
in the literature). The answers for some of these definite integrals involving Jacobi elliptic functions are 
given in Appendix G. 

Actually, we can even simplify several indefinite integrals. In fact, by starting from any local identity 
we can obtain the indefinite integral of its left hand side in terms of the well known integrals [31 E] of 
sn n (x), dn n (x), cn n (x). The only exceptions are those MI-II local identities in which the right hand side 
has a term proportional either to [Z(x + a) + Z(x — a) — 2Z(x)] or to [Z(x + a) — Z(x — a)]. We show below 
that in that case the indefinite integral of the left hand side also has terms containing indefinite elliptic 
integrals of the first, second and third kind [31 E]. 

We start from the local identity 

dn 2 (x)dn(x + a) = Bdn(x + a) + ds(a)ns(a)dn(x) — mcs(a)sn(x)cn(x) , B = — cs 2 (a) , (59) 

which is obtained by adding the identities 1)26(1 and (|32j) . On integrating both sides with respect to x and 
using the known integral of dn(x) [5] we then obtain 

/ dn^)dn(, + a) dx = £am ( * + a) + ds(a)ns(a)am ( ,) + cs ( a)dn(x) , (60) 
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where / dn(x) dx = am(x) = sin (sn(x)) = iln [cn(x) — isn(x)]. On multiplying both sides of eq. by 
dn 2 (x) and using eq. (|60|) we then find that 

+ds(a)ns(a 



1 + 

2 



, s r, r s -i / \ cs(a) , o. . mdsfa)ns(a) . . , . 

am(x) + 5cs(a)dn(x) H — dn d (x) H — — — sn(x)cn(x) . (61) 
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Proceeding in this way, one can find the indefinite integral of dn 2n (x)dn(x + a). One can also derive a 
recursion relation relating the various integrals. In particular, using eq. Q35|) one can show that 

dn 2n (x)dn(x + a) = J Bdn 2n_2 (x)dn(2; + a) + ds(a)ns(a)dn 2n_1 (x) - mcs(a)sn(x)cn(x)dn 2n ~ 2 (x) . (62) 

On integrating both sides with respect to x, we find a recursion relation relating various integrals: 

/„ = SJ„_i + CS ( Q ) dn 2 "-V) + ds(a)ns(a) / dn 2 "" 1 ^) dx ; I k = [ dn 2fc (x)dn(x + a) dx. (63) 
(2n — 1) J J 

Note that the integral of any power of dn(x) (as well as sn(x), cn(x)) is known in principle We might 
also add that once the integral of say dn 2n (x)dn(x + a) is known then the integral of dn(x)dn 2n (x + a) is 
obtained from it by simply replacing x by x — a followed by a — > —a. 

Using the above procedure one can obtain indefinite integrals of the left hand sides of all the local 
identities given in this paper except for those MI-II local identities where the combination [Z(x+a)— Z(x— a)] 
or [Z(x + a) + Z(x — a) — 2Z(x)] also occurs on the right hand side. Let us now explain how to handle these 
integrals. 

We start from the local identity (JT7|) . On using eq. © it is easily shown that 

, / sr, / s , / m 2dn(a)[l — msn 2 (x)l , ,. 

dn(x)[dn(x + a) + dn(x - a)] = - V Jl Tl - [ ] . (64) 

1 — msn z (a)sn z (x) 

On integrating both sides of eq. (|17|) over x, using eq. (|64|) and well known integrals (see integrals 336.01 
and 337.01 of (0)) we finally find that 

J dn(x)[dn(x + a) + dn(x — a)] dx = 2[dn(a) — cs(a)Z(a)]x + cs(a) J [Z(x + a) — Z(x — a)] dx 

= 2ds(a)ns(a)F(am x, k) — 2dn(a)cs 2 (a)n(am x, /c 2 sn 2 (a), k) , (65) 
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where k 2 = m . Here F(am x, k) and Il(am x, & 2 sn 2 (a), k) are indefinite elliptic integrals of first and third 
kind respectively. Similarly, on using eqs. © and (|18j) it is easy to show that 

fM*)[M*+*)~M*-a)] d, = 08(a) /[Z(s + a) + Z(s-a)-2Z(s)] dx = cs ( a) ln[l-msn 2 (a)sn 2 (*)], (66) 
and hence 



dn(x)dn(x + a) dx = [dn(a) — cs(a)Z(a)]x — cs(a) / [Z(x + a) — Z(x)] dx 
= ds(a)ns(a)-F(am x, k) — dn(a)cs 2 (a)IT(am x, k 2 sn 2 (a), k) + (l/2)cs(a) ln[l — msn 2 (a)sn 2 (x)]. (67) 

We now show that using eq. (|67|) one can obtain the indefinite integral of the left hand side of any 
local MI-II identity. As an illustration, consider the local identity (|2U8|) . It is easy to show from here that 

dn 2n (x)dn 2 (x + a) = £dn 2n " 2 (x)[dn 2 (x) + dn 2 (x + a)] - (1-m )dn 2n " 2 (x) + 25 n " 1 ds(a)ns(a)dn(x + a)dn(x) 

n-1 

+2ds(a)ns(a ) [ds(a)ns(a)dn(x) - m cs(a)sn(2;)cn(x)] ^ B^ 1 [dn^)] 2 ^"^" 1 , (68) 

k=l 

where B = — cs 2 (a) and n > 2. On integrating both sides of this equation, we get the recursion relation 
I n = B^ + J Bdn 2n (x) dx-{l-m) J dn 2n ~ 2 (x) dx + 2£ n - 1 ds(a)ns(a) J dn(x + a)dn(x) dx 

n-1 n 2n-2k( \ n-1 

+cs(a)ds(a)ns(a) ^ B k ~ x — + 2ds 2 (a)ns 2 (a) ^ / B^dn^-^ix) dx , (69) 

k=i n ~ k k=i J 

where I\ is easily obtained by using the integral (|BTj) and the local identity (|166j) . We find 
I^Jdn 2 (x)dn\x + a)dx=-c S 2 ^ 

(70) 

where E(am x, k) is the indefinite elliptic integral of the second kind. 

Recursion relations for several indefinite elliptic integrals are given in Appendix H, where the right 
hand side is in terms of Jacobi elliptic functions, their integrals, and indefinite elliptic integrals of the first, 
second and third kind. 

6 Continuum Limit of Local and Cyclic Identities 

We now study what happens to the local identities as a — > 0. Although the identities are not valid at a = 0, 
this limit leads to well known nonlinear ordinary differential equations satisfied by the Jacobian elliptic 
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functions. Thus we will see that the local identities may be viewed as exact discretization of these differential 
equations. This provides our justification for calling these identities "local". This is to be contrasted with 
the cyclic identities that are exact discretization of integral identities. Just as the differential equation can 
be integrated, the local identities are simply summed to produce the cyclic identities. 
Take the simple type I identity: 

dn 2 (z)[dn(z + a) + dn(z - a)] = 2ds(a)ns(a)dn(z) - cs 2 (a)[dn(z + a) + dn(z - a)] . (71) 

Since cs 2 (a) has a pole of order two at a = 0, we expand to second order in a: 

d 2 

dn(z + a) + dn(z - a) = 2dn(z) + a 2 -^dn(z) + 0(a 3 ) . (72) 

Using 

lim [ds(a)ns(a) — cs 2 (a)] = 1 — — , (73) 
then leads to the limiting differential equation: 

(2-m)y-0 = 2y 3 , (74) 

with y = dn(z) . Of the many applications of such differential equations, the most straightforward one 
would be the interpretation of this as Newton's equation of motion for a particle in a one-dimensional 
double well potential, with z taking the role of time. 

It is well known that finite difference versions of such nonlinear one-dimensional problems tend to 
exhibit chaos, and therefore the difference equation does not have analytical solutions in terms of the 
Jacobi elliptic functions. However we precisely achieve this when converting this differential equation into 
a finite difference equation in the following manner: 

(2 _ m)y(z) b(, + A) + »(,-A)- W ] = y2(zMz + A) + y(z A)] (?5) 

Of course depending on the smallness of A, the above discrete equation will only be an approximation to 
the actual solution dn(z). However, replacing 2 — m in the difference scheme with 2[ds(A)ns(A) — cs 2 (A)] 
and the 1/A 2 factor multiplying the second difference with cs 2 (A) leads to an exact difference scheme: 

2[ds(A)ns(A) - cs 2 (A)]?/ n - cs 2 (A)[y n+1 + y n _i - 2y n ] = y 2 [y n +i + Vn-i] , (76) 

20 



where y(n) = y{z + nA), which is identical to the local identity in eq. Q71JI when we put A = a. This is of 
course just the reverse of the small a limit. 

The local identity implies that y n = dn(a n) is a solution of this difference equation with the initial 
condition yo = 1, y\ = dn(a). The modulus parameter is m in all the elliptic functions involved. This is 
analogous to the differential equation eq. (|74|) whose solution is y(z) = dn(z), when the initial conditions 
are y(0) = l,y'(0) = 0, where the prime denotes differentiation. It may be noted that the local identity 
in eq. (|43|) also has the same continuum limit as the above, and that the difference equation is therefore 
genuinely second order rather than first. All of the two-point local identities limit to well known differential 
equations, or to those that can be derived from these. 

The cyclic identities limit to definite integrals, but we may, with little effort, also generalize to indefinite 
integrals. By way of illustration again consider the type I local identity in eq. (|71|) . Let Z{ = zq + A (i — l)/p, 
with i = 1, . . . , p. Here A is arbitrary and p is an integer, Zi are the sample points between zq and z = zq+A. 

Chaining the local identities at these sample points together we get: 

v p 
^2dn 2 (zi)[dn(z i+1 ) + dn(zj_i)] = [2ds(a)ns(a) - cs 2 (a)] dn(zj) 
i=i i=i 

-cs 2 (a)[dn(2;p + i) - dn(zi) + dn(z ) - dn(z p )] , (77) 

with a = A I p. This becomes a cyclic identity when A = 2K in which case the "end correction" that is the 
last term in the above equation vanishes due to the periodicity of dn(z). Multiply both sides by A/p and 
then take the a — > (or equivalently p — > oo) limit. The end corrections are rewritten and evaluated as 
follows: 

A p 

cs 2 (a)[dn(z p+ i) - dn(zi) + dn(z ) - dn(z p )\ -> — -[dn(z p+ i) - dn(z p ) + dn(z ) - dn(zi)] (78) 

p A 

— > m sn(z)cn(z) — m sn(zo)cn(zo) • 

Here we have used that the derivative of dn(z) is -msn(z)cn(z). Thus the entire identity in eq. (|77|) limits 
effectively to the indefinite integral: 

J dn 3 (z) dz = -(2 - m) J dn(z) dz + ysn(z)cn(z) , (79) 

which is a standard identity, for instance eq. (314.03) of Byrd and Friedman [Bj. Thus the cyclic identities 
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are exact Riemann sums of integral identities while local identities are exact discretization of differential 
equations. 

In the case of three-point local identities (those which have the elliptic functions evaluated at three 
points, z, z + a and z + a') there are more than one limiting cases. The case a' — * a or a' —* leads to 
two-point local identities involving z and z + a. For instance starting from the local identity in eq. (|93j) 
and taking the a — > a' limit or the a' — ► limit leads essentially to the local two-point identity in eq. (|43|). 

7 Comments and Discussion 

In this paper, we have proved a wide class of local identities, using which we have obtained corresponding 
cyclic identities with arbitrary weights. These results are easily extended in several directions. For example, 
we could evaluate these identities at points separated by gaps of T/p with imaginary or complex period 
T thereby obtaining corresponding local identities for pure imaginary as well as complex shifts. Secondly, 
we can also obtain similar identities for all 9 auxiliary Jacobi functions like ns(x,m) as well as 6 ratios of 
Jacobi functions like cn(x)dn(x)/sn(x). Further, by following the procedure in ref. II, we can readily write 
down the corresponding local as well as cyclic identities for Weierstrass functions, as well as for the ratio 
of any two of the four Jacobi theta functions. 

Many of these identities have nontrivial m = 0, 1 limits. For example, a cyclic identity derived in ref. 

II, which is valid for any odd integer p, and I < p, is 

p r (»-i)/2 (Z-i)/2 i i p 

Y J d j d j+r ...d m ^y= J] ^N + ^-ljC- 1 )/ 2 ^ J] cs([n-k)a) J> , (80) 

k=l k=l n=l,nj^k 



3=1 

while for I = p, we have the simpler identity 



3=1 



{p - 1)/2 ,2Kn p 



Y[ dj = n c S 2 (— ^d,, (si) 

3=1 ™=1 3=1 

where a = rlKjp. At m =0, these identities reduce to interesting trigonometric identities 

(Z-l)/2 (J-l)/2 I 

1= J] cot 2 (™) + 2(-l)('- 1 )/ 2 2 I] ootan-*]^), (82) 

h=\ P k=l n=l,n^k " 

1 (P ^ /2 ,,n, 



p , p 

r n=l ' 
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Actually, identity is also valid for even p and r = 1, provided I < {jp + 2)/2. 

For the special case of / = 3, one can write down cyclic identities for products of three dn's at arbitrary 
separation (in units of 2K/p) for both even and odd p. In particular, using the local identity (|93j) . one can 
immediately write down cyclic identities for combinations like Yjj=\ djdj+ r dj+ s with r and s being unequal 
but arbitrary otherwise. In the limit m = 0, summing of all such independent cyclic identities yields the 
following remarkable trigonometric identities: 

(p-l)(p-2) = g cot2(g) (84) 

Similarly, many of the other local identities we have derived in this paper also reduce to interesting 
trigonometric identities in the limit m = 0. 

Finally, note that even though our local identities have been derived assuming that a is an arbitrary 
constant, the identities are also valid when a is any function of x. For example, on calling x + a = b, the 
local identities (|§U)) and (|§T|) give generalized addition theorems 

dn(a — 6)sn(a)sn(6) + cn(a)cn(6) = cn(a — b) , (85) 

dn(a — 6)sn(a)cn(6) — cn(a)sn(6) = — dn(a)sn(a — b) , (86) 

which in the limit m = reduce to the well known addition theorems for the trigonometric functions. 
AK and US gratefully acknowledge support from the U.S. Department of Energy. 
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Appendix A: Local Identities of Rank 2 

Rank 2 identities with 2 distinct arguments (x,x + a). 



dn(x)dn(x + a) = dn(a) + cs(a) [Z(x + a) - Z(x) - Z(a)] , (87) 

msn(x)sn(x + a) = -ns(a) [Z(x + a) - Z(x) - Z(a)] , (88) 

mcn(x)cn(x + a) = mcn(a) + ds(a) [Z(x + a) — Z(x) — Z(a)] , (89) 

dn(x)sn(x + a) = ns(a)cn(x) — cs(a)cn(x + a) , (90) 

dn(x)cn(x + a) = — ds(a)sn(x) + cs(a)sn(x + a) , (91) 

msn(x)cn(x + a) = ds(a)dn(x) — ns(a)dn(x + a) . (92) 

Appendix B: Local Identities of Rank 3 

Rank 3 identities with 3 distinct arguments (x, x + a, x + a'). 

dn(x)dn(x + a)dn(x + a') = — cs(a)cs(a')dn(x) — cs(a)cs(a — a')dn(x + a) + cs(a')cs(a — a')dn(x + a') , (93) 

msn(x)sn(x + a)sn(x + a') = ns(a)ns(a')sn(x) + ns(a)ns(a — a')sn(x + a) — ns(a')ns(a — a')sn(x + a') , (94) 

mcn(x)cn(x+a)cn(x+a') = — ds(a)ds(a')cn(x) — ds(a)ds(a — a')cn(x + a)+ds(a')ds(a — a')cn(x+a') , (95) 

dn(x)dn(x + a)sn(x + a') = — cs(a)ns(a')sn(x) — cs(a)ns(a — a')sn(x + a) + cs(a')cs(a — a')sn(x + a') , (96) 

dn(x)dn(x + a)cn(x + a') = — cs(a)ds(a')cn(x) — cs(a)ds(a — o')cn(x + a) + cs(a')cs(a — a')cn(x + a') , (97) 

msn(x)sn(x + a)dn(x + a') = ns(a)cs(a')dn(x) + ns(a)cs(a — a')dn(x + a) — ns(a')ns(a — a')dn(x + a') , (98) 

msn(x)sn(x + a)cn(x + a') = ns(a)ds(a')cn(x) + ns(a)ds(a — a')cn(x + a) — ns(a')ns(a — a')cn(x + a') , (99) 

mcn(x)cn(x+a)dn(x+a') = — ds(a)cs(a')dn(x)— ds(a)cs(a— a')dn(x+a)+ds(a')ds(a— a')dn(x+a') , (100) 

mcn(x)cn(x+a)sn(x+a') = — ds(a)ns(a')sn(x) — ds(a)ns(a— a')sn(x+a)+ds(a / )ds(a— a')sn(x+a') , (101) 
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mdn(x)sn(x + a)cn(x + a ) = — ds(a — a'){dn(a) + cs(a) [Z(x + a) — Z(x) — Z(a)]} 

+ns(a - o'){dn(a') + cs(a') [Z(x + a') - Z(x) - Z(a')] } • (102) 



Rank 3 identities with 2 distinct arguments (x,x + a). 

As explained in the text, local identities for say m cn(x)sn(x + a)dn(x + a) and m sn(x)dn(x)cn(x + a) 
are related to each other by x — > x — a followed by a — > — a and hence only one of these identities is given 
below. 



dn 2 (x)dn(x + a) = — cs 2 (a)dn(x + a) + ds(a)ns(a)dn(x) — mcs(a)cn(x)sn(x) , (103) 

msn 2 (x)sn(x + a) = ns 2 (a)sn(x + a) — cs(a)ds(a)sn(x) — ns(a)cn(x)dn(x) , (104) 

mcn 2 (x)cn(x + a) = — ds 2 (a)cn(x + a) + cs(a)ns(a)cn(x) — ds(a)sn(x)dn(x) , (105) 

dn(x)sn(x)dn(x + a) = — cs(a)ns(a)sn(x + a) + ds(a)ns(a)sn(x) + cs(a)cn(x)dn(x) , (106) 

dn(x)cn(x)dn(x + a) = — cs(a)ds(a)cn(x + a) + ds(a)ns(a)cn(x) — cs(a)sn(x)dn(x) , (107) 

mdn(x)sn(x)sn(x + a) = cs(a)ns(a)dn(x + a) — cs(a)ds(a)dn(x) + mns(a)cn(x)sn(x) , (108) 

msn(x)cn(x)sn(x + a) = ds(a)ns(a)cn(x + a) — cs(a)ds(a)cn(x) + ns(a)sn(x)dn(x) , (109) 

mdn(x)cn(x)cn(x + a) = — cs(a)ds(a)dn(x + a) + cs(a)ns(a)dn(x) — mds(a)cn(x)sn(x) , (HO) 

msn(x)cn(x)cn(x + a) = — ds(a)ns(a)sn(x + a) + cs(a)ns(a)sn(x) + ds(a)cn(x)dn(x) , (HI) 

mdn(x)sn(x)cn(x + a) = — ds(a) — cs(a)ns(a) [Z(x + a) — Z(x) — Z(a)] + ds(a)dn 2 (x) , (H2) 
mcn(x)dn(x)sn(x + a) = — ds(a)dn(a) — cs(a)ds(a) [Z(x + a) — Z(x) — Z(a)] + ns(a)dn 2 (x) , (H3) 

msn(x)cn(x)dn(x + a) = — cs(a) — ds(a)ns(a) [Z(x + a) — Z(x) — Z(a)] + cs(a)dn 2 (x) . (H4) 



Appendix C: Local Identities of Rank 4 
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Rank 4 identities with 4 distinct arguments (x, x + a, x + a', x + a"). 

dn(x)dn(x + a)dn(x + a')dn(x + a") = -cs(a)cs(a'){dn(a") + cs(o")[Z(x + a") - Z(x) - Z(a")]} 
-cs(a)cs(a - a'){dn(a" - a) + cs(a" - a)[Z(x + a") - Z(x + a) - Z(o" - a)]} 
+cs(a')cs(a - a'){dn(a" - a') + cs(a" - a')[Z(x + a") - Z(x + a') - Z(a" - a')]} , (115) 

m 2 sn(x)sn(x + a)sn(x + a')sn(x + a") = — ns(a)ns(a')ns(a")[Z(x + a") — Z(x) — Z(a")] 



m 2 cn(a;)cn(x + a)cn(x + a')cn(x + a") = -ds(o)ds(a'){mcn(a") + ds(a")[Z(x + a") - Z(x) - Z(a")]} 
-ds(a)ds(a - a'){mcn(o" - a) + ds(a" - a)[Z(x + a") - Z(x + a) - Z(a" - a)]} 
+ds(a')ds(a - a'){mcn(a" - a') + ds(a" - a')[Z(x + a") - Z(x + a') - Z(a" - a')]} , (117) 

dn(x)dn(x + a)dn(x + a')sn(x + a") = — cs(a)cs(a / )[ns(a // )cn(x) — cs(a")cn(x + a")] 



-ns(a)ns(a - a')ns(a" - a)[Z(x + a") - Z(x + a) - Z(o" - a)] 



+ns(a')ns(a — a')ns(a" — a')[Z(x + a") — Z(x + a') — Z(a" — a')] , 



(116) 



— cs(a)cs(a — a')[ns(a" — a)cn(x + a) — cs(a" — a)cn(x + a")] 



+cs(a')cs(a — a')[ns(a" — a')cn(x + a') — cs(a" — a')cn(x + a")] , 



(118) 



dn(x)dn(x + a)dn(x + a')cn(x + a") = cs(a)cs(a')[ds(a")sn(x) — cs(a")sn(x + a")] 



+cs(a)cs(a — a')[ds(a" — a)sn(x + a) — cs(a" — a)sn(x + a")] 



— cs(a')cs(a — a')[ds(a" — a')sn(x + a') — cs(a" — a')sn(x + a")] , 



(119) 



msn(x)sn(x + a)sn(x + a')dn(x + a") = ns(a)ns(a / )[cs(a")cn(x) — ns(a")cn(x + a")] 



+ns(a)ns(a — a')[cs(a" — a)cn(x + a) — ns(a" — a)cn(x + a")] 



— ns(a')ns(a — a')[cs(a" — a')cn(x + a') — ns(a" — a')cn(x + a")] , 



(120) 



m 2 sn(x)sn(x + a)sn(x + a')cn(x + a") = ns(a)ns(a')[ds(a")dn(x) — ns(a")dn(x + a")] 



+ns(a)ns(a — a')[ds(a" — a)dn(x + a) — ns(a" — a)dn(x + a")] 



-ns(a')ns(a - o')[ds(a" - a')dn(x + a') - ns(a" - a')dn(x + a")] 



(121) 
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mcn(x)cn(x + a)cn(x + a')dn(x + a") = ds(a)ds(a')[cs(a")sn(x) — ds(a")sn(x + a")] 



+ds(a)ds(a — a')[cs(a" — a)sn(x + a) — ds(a" — a)sn(x + a")} 



-ds(a')ds(a - a')[cs(a" - a')sn(x + a') - ds(a" - a')sn(x + a")] 



(122) 



m 2 cn(x)cn(x + a)cn(x + a')sn(x + a") = — ds(a)ds(a')[ns(a")dn(x) — ds(a")dn(x + a")] 



— ds(a)ds(a — a')[ns(a" — a)dn(x + a) — ds(a" — a)dn(x + a")] 



+ds(a')ds(a - o')[ns(a" - a')dn(x + a') - ds(a" - o')dn(a; + a")} 



(123) 



msn(i)dn(x + a)sn(x + a')dn(x + a") = cs(a)ns(a / ){dn(a") + cs(a")[Z(x + a") - Z(x) - Z(a")]} 
+ns(a)ns(a - a'){dn(a" - a) + cs(a" - a)[Z(x + a") - Z(x + a) - Z(a" - a)]} 
-ns(a')cs(a - a ){dn(a" - a') + cs(a" - a')[Z(x + a") - Z(x + a') - Z(a" - a')]} , (124) 

mcn(x)dn(x + a)cn(x + a')dn(x + a") = -cs(a)ds(a'){dn(a") + cs(a")[Z(x + a") - Z(x) - Z(a")]} 
-ds(a)ds(a - a'){dn(a" - a) + cs(a" - a)[Z(x + a") - Z(x + a) - Z(o" - a)]} 
+ds(a')cs(a - a'){dn(a" - a) + cs(a" - a')[Z(x + a") - Z(x + a) - Z(a" - a')]} , (125) 

m 2 sn(x)cn(x + a)cn(x + a')sn(x + a") = ds(a)ds(a')ns(a")[Z(x + a ") _ z(x) - Z(a")] 



+ns(a)ds(a - a')ns(o" - a)[Z(x + a") - Z(x + a) - Z(a" - a)] 



-ns(a')ds(o - o')ns(o" - a')[Z(x + a") - Z(x + a') - Z(a" - a')] 



(126) 



mcn(x)dn(x + a)dn(x + a')sn(x + a") = — cs(a)cs(a')[ns(a")dn(x) — ds(a")dn(x + a")] 



— ds(a)cs(a — a') [ns(a — a)dn(x + a) — ds(a — a)dn(x + a )] 



+ds(a')cs(a - a')[ns(a" - a')dn(x + a') - ds(a" - a')dn(x + a")] , 



(127) 



msn(x)dn(x + a)sn(x + a')cn(x + a") = — cs(a)ns(a')[ds(a")sn(x) — cs(a")sn(x + a")] 



— ns(a)ns(a — a')[ds(a" — a)sn(x + a) — cs(a" — a)sn(x + a")] 



+ns(a')cs(a — a')[ds(a" — a')sn(x + a') — cs(a" — a')sn(x + a")] 



(128) 
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mcn(x)dn(x + a)cn(x + a')sn(x + a") = — cs(a)ds(a')[ns(a")cn(x) — cs(a // )cn(x + a")] 
— ds(a)ds(a — a')[ns(a" — a)cn(x + a) — cs(a" — a)cn(x + a")} 
+ds(a')cs(a - a') [ns(a" - a')cn(x + a') - cs(a" - a')cn(x + a")} . (129) 



Rank 4 identities with 3 distinct arguments (x,x + a,x + a'). 



dn 2 (x)dn(x + a)dn(x + a') = — cs(a)cs(a — a'){dn(a) + cs(a)[Z(x + a) — Z(x) — Z(a)]} 

+cs(a')cs(a — a'){dn(a') + cs(a')[Z(x + a') — Z(x) — Z(a')]} — cs(a)cs(a ')dn 2 (x) (130) 

m 2 sn 2 (x)sn(x + a)sn(x + a') = — ns 2 (a)ns(a — a')[Z(x + a) — Z(x) — Z(a)] 

+ns 2 (a')ns(a - a')[Z(x + a') - Z(x) - Z(a')] + mns(a)ns(a')sn 2 (x) , (131) 

m 2 cn 2 (x)cn(x + a)cn(x + a') = — ds(a)ds(a — o'){mcn(a) + ds(a)[Z(x + a) — Z(x) — Z(a)]} 

+ds(a')ds(a - a / ){mcn(a / ) + ds(a')[Z(x + a') - Z(x) - Z(a')]} - mds(a)ds(a')cn 2 (x) ,(132) 

mdn(x)sn(x)dn(x + a)sn(x + a') = cs(a)ns(a)ns(a — a')[Z(x + a) — Z(x) — Z(a)] 

-cs(a )ns(o')cs(o - a')[Z(x + a) - Z(x) - Z(a )] - mcs(a)ns(a')sn 2 (x) , (133) 

mdn(x)cn(x)dn(x + a)cn(x + a') = — cs(a)ds(a — a'){mcn(a) + ds(a)[Z(x + a) — Z(x) — Z(a)]} 

+cs(a')cs(a — a '){mcn(a') + ds(a')[Z(x + a) - Z(x) - Z(a )]} - mcs(a)ds(a')cn 2 (x) , (134) 

m 2 sn(x)cn(x)sn(x + a)cn(x + a') = ds(a)ns(a)ds(a — a')[Z(x + a) — Z(x) — Z(a)] 

— ds(a')ns(a')ns(a — a ')[Z(x + a ) - Z(x) - Z(a )] - mns(a)ds(a')sn 2 (x) , (135) 



dn(x) 2 dn(x + a)sn(x + a) = — cs(a)ns(a — a')[ns(a)cn(x) — cs(a)cn(x + a)] 

+cs(a')cs(a — a')[ns(a')cn(x) — cs(a')cn(x + a')] — cs(a)ns(a')sn(x)dn(x) , 
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(136) 



dn(x) 2 dn(x + a)cn(x + a') = cs(a)ds(a — a')[ds(a)sn(x) — cs(a)sn(x + a)] 

— cs(a')cs(a — a')[ds(a')sn(x) — cs(a')sn(x + a')] — cs(a)ds(a')cn(x)dn(x) , 



(137) 



m 2 cn(x) 2 sn(x + a)cn(x + a') = — ds(a)ds(a — a')[ns(a)dn(x) — ds(a)dn(x + a)] 

+ds(a')ns(a — a')[ns(a / )dn(x) — ds(a')dn(x + a')] — mns(a)ds(a / )cn(x)sn(x) , (138) 



mdn(x)sn(x)sn(x + a)cn(x + a) = — ns(a)ds(a — a')[ds(a)sn(x) — cs(a)sn(x + a)] 

+ns(a / )ns(a — a')[ds(a')sn(x) — cs(a')sn(x + a')] + ns(a)ds(a')cn(x)dn(x) , (139) 



mdn(x)sn(x)dn(x + a)cn(x + a') = — cs(a)ds(a — a')[ds(a)dn(x) — ns(a)dn(x + a)] 

+cs(a')cs(a — a / )[ds(a')dn(x) — ns(a')dn(x + a')] — mcs(a)ds(a')cn(x)sn(x) , (140) 



dn(x)sn(x)dn(x + a)dn(x + a') = — cs(a)cs(a — a')[cs(a)cn(x) — ns(a)cn(x + a)] 

+cs(a')cs(a — a')[cs(a')cn(x) — ns(a')cn(x + a')} — cs(a)cs(a')sn(x)dn(x) , (141) 



mdn(x)sn(x)sn(x + a)sn(x + a') = ns(a)ns(a — a')[ns(a)cn(x) — cs(a)cn(x + a)] 

— ns(a')ns(a — a / )[ns(a')cn(x) — cs(a')cn(x + a')] + ns(a)ns(a')sn(x)dn(x) , (142) 



mdn(x)sn(x)cn(x + a)cn(x + a) = — ns(a)ds(a — a')[ns(a)cn(x) — cs(a)cn(x + a)] 
+ns(a')ds(a — a')[ns(a')cn(x) — cs(a')cn(x + a')] — ds(a)ds(a / )sn(x)dn(x) , 



(143) 



mdn(x)cn(x)sn(x + a)cn(x + a) = — ds(a)ds(a — a')[ns(a)cn(x) — cs(a)cn(x + a)] 
+ds(a')ns(a — a / )[ns(a')cn(x) — cs(a')cn(x + a')] — ns(a)ds(a')sn(x)dn(x) , 



(144) 



mdn(x)cn(x)dn(x + a)sn(x + a') = — cs(a)ns(a — a')[ns(a)dn(x) — ds(a)dn(x + a)] 

+cs(a')cs(a — a')[ns(a')dn(x) — ds(a')dn(x + a')] — mcs(a)ns(a')cn(x)sn(x) , (145) 
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dn(x)cn(x)dn(x + a)dn(x + a) = cs(a)cs(a — a')[cs(a)sn(x) — ds(a)sn(x + a)] 

— cs(a')cs(a — a')[cs(a / )sn(x) — ds(a')sn(x + a')] — cs(a)cs(a')cn(x)dn(x) , (146) 

mdn(x)cn(x)sn(x + a)sn(x + a) = — ds(a)ns(a — a')[ds(a)sn(x) — cs(a)sn(x + a)] 

+ds(a')ns(a — a')[ds(o')sn(x) — cs(a')sn(x + a')] + ns(a)ns(a')cn(x)dn(x) , (147) 

mdn(x)cn(x)cn(x + a)cn(x + a') = ds(a)ds(a — a')[ds(a)sn(x) — cs(a)sn(x + a)] 

— ds(a')ds(a — a')[ds(a')sn(x) — cs(a')sn(x + a')] — ds(a)ds(a')cn(x)dn(x) , (148) 

m 2 sn(x)cn(x)sn(x + a)sn(x + a') = ds(a)ns(a — a')[ds(a)dn(x) — ns(a)dn(x + a)] 

— ds(a')ns(a — a')[ds(a')dn(x) — ns(a')dn(x + a')] + mns(a)ns(a')sn(x)cn(x) , (149) 

msn(x)cn(x)dn(x + a)dn(x + a') = — ns(a)cs(a — a')[ns(a)dn(x) — ds(a)dn(x + a)] 

+ns(a')cs(a — a')[ns(a')dn(x) — ds(a')dn(x + a')] — mcs(a)cs(a')cn(x)sn(x) , (150) 

m 2 sn(x)cn(x)cn(x + a)cn(x + a') = — ns(a)ds(a — a')[ns(a)dn(x) — ds(a)dn(x + a)] 

+ns(a')ds(a — a')[ns(a')dn(x) — ds(a')dn(x + a')] — mds(a)ds(a')cn(x)sn(x) , (151) 

msn(x)cn(x)dn(x + a)sn(x + a') = — ns(a)ns(a — a')[cs(a)sn(x) — ds(a)sn(x + a)] 

+ns(a')cs(a — a')[cs(a')sn(x) — ds(a')sn(x + a')} + cs(a)ns(a')cn(x)dn(x) , (152) 

msn(x)cn(x)dn(x + a)cn(x + a') = — ds(a)ds(a — a')[cs(a)cn(x) — ns(a)cn(x + a)] 

+ds(a')cs(a — a')[cs(a / )cn(x) — ns(a')cn(x + a')] — cs(a)ds(a / )sn(x)dn(x) . (153) 

Rank 4 identities with 2 distinct arguments (x, x + a). 

mdn(x)sn(x)cn(x)dn(x + a) = — cs(a)[ds 2 (a) + l]dn(x) + cs(a)ds(a)ns(a)dn(x + a) 

+ cs(a)dn 3 (x) + mds(a)ns(a)cn(x)sn(x) , (154) 
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mdn(x)sn(x)cn(x)sn(x + a) = — ns(a)[ds 2 (a) — l]sn(x) + cs(a)ds(a)ns(a)sn(x + a) 

— mns(a)sn 3 (x) — cs(a)ds(a)cn(x)dn(x) , (155) 

mdn(x)sn(x)cn(x)cn(x + a) = — ds(a)[m + cs 2 (a)]cn(x) + cs(a)ds(a)ns(a)cn(x + a) 

+mds(a)cn 3 (x) + cs(a)ns(a)sn(x)dn(x) , (156) 

dn 3 (x)dn(x + a) = — mcs(a)sn(x)cn(x)dn(x) + ds(a)ns(a)dn 2 (x) 

-cs 2 (a)dn(a)-cs 3 (a)[Z(x + a)-Z(x)-Z(a)] , (157) 

dn 3 (x)sn(x + a) = cs 3 (a)cn(x + a) — ns(a) ^2cs 2 (a) — ds 2 (a)) cn(x) 

+ cs(a)ds(a)sn(x)dn(x) + mns(a)cn 3 (x) , (158) 

dn 3 (x)cn(x + a) = — cs 3 (a)sn(x + a) — ds(a) ^2 — ns 2 (a)^ sn(x) 

+ cs(a)ns(a)cn(x)dn(x) + mds(a)sn 3 (x) , (159) 

msn 3 (x)dn(x + a) = — ns 3 (a)cn(x + a) + + ns 2 (a)^ cs(a)cn(x) 

— ds(a)ns(a)sn(x)dn(x) — mcs(a)cn 3 (x) , (160) 

m 2 sn 3 (x)sn(x + a) = — mns(a)sn(x)cn(x)dn(x) + cs(a)ds(a)dn 2 (x) 

-cs(a)ds(a)-ns 3 (a)[Z(x + a)-Z(x)-Z(a)] , (161) 

m 2 sn 3 (x)cn(x + a) = ^1 + ns 2 (a)^ ds(a)dn(x) — ns 3 (a)dn(x + a) 

— mcs(a)ns(a)sn(x)cn(x) — ds(a)dn 3 (x) , (162) 

mcn 3 (x)dn(x + a) = — ds 3 (a)sn(x + a) + ^2ds 2 (a) — ns 2 (a)) cs(a)sn(x) 

+ ds(a)ns(a)cn(x)dn(x) + mcs(a)sn 3 (x) , (163) 
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m 2 cn 3 (x)sn(x + a) = ds 3 (a)dn(x + a) — ^2ds 2 (a) — cs 2 (a)) ns(a)dn(x) 
+ mcs(a)ds(a)sn(x)cn(x) + ns(a)dn 3 (x) , 



(164) 



m 2 cn 3 (x)cn(x + a) = — mds(a)sn(x)cn(x)dn(x) + cs(a)ns(a)dn 2 (x) 

-cs(o)ns(a)[l - m + mdn 2 (a)] - ds 3 (a) [Z(x + a) - Z(x) - Z(a)] , (165) 

dn 2 (x)dn 2 (x + a) = -cs 2 (a)[dn 2 (x) + dn 2 (x + a)] + [ds 2 (a) + cs 2 (a)] 

+2cs(a)ds(a)ns(a) [Z(x + a) - Z(x) - Z(a)] , (166) 

mdn 2 (x)sn(x + a)cn(x + a) = cs(a)[ds 2 (a) + ns 2 (a)]dn(x) — 2cs(a)ds(a)ns(a)dn(x + a) 

— mcs 2 (a)sn(x + a)cn(x + a) — mds(a)ns(a)cn(x)sn(x) , (167) 

msn 2 (x)dn(x + a)cn(x + a) = ns(a)[cs 2 (a) + ds 2 (a)]sn(x) — 2cs(a)ds(a)ns(a)sn(x + a) 

+ ns 2 (a)cn(x + a)dn(x + a) + cs(a)ds(a)cn(x)dn(x) , (168) 

mcn 2 (x)dn(x + a)sn(x + a) = ds(a)[cs 2 (o) + ns 2 (a)]cn(x) — 2cs(a)ds(a)ns(a)cn(x + a) 

— ds 2 (a)sn(x + a)dn(x + a) — cs(a)ns(a)sn(x)dn(x) , (169) 



m 2 sn(x)cn(x)sn(x + a)cn(x + a) = ds(a)ns(a) dn 2 (x) + dn 2 (x + a) 

- ds 2 (a) + ns 2 (a)l {dn(a) + cs(a) [Z(x + a) - Z(x) - Z(a)]} , (170) 

mdn(x)sn(x)dn(x + o)sn(x + a) = — cs(a)ns(a)[l + dn 2 (a)] + cs(a)ns(a) dn 2 (x) + dn 2 (x + a) 

-ds(a)[cs 2 (a) + ns 2 (a)][Z(x + a)-Z(x)-Z(a)] , (171) 



mdn(x)cn(x)dn(x + a)cn(x + a) = 2cs(a)ds(a) — cs(a)ds(a) dn 2 (x) + dn 2 (x + a) 

+ns(a)[ds 2 (a) + cs 2 (a)] [Z(x + a) - Z(x) - Z(a)] , (172) 
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mdn(x)cn(x)dn(x + a)sn(x + a) = ds(a)[cs 2 (a) + ns 2 (a)]dn(x) — ns(a)[cs 2 (a) + ds 2 (a)]dn(x + a) 

— mcs(a)ds(a)sn(x + a)cn(x + a) — mcs(a)ns(a)cn(x)sn(x) , (173) 

mdn(x)cn(x)sn(x + a)cn(x + a) = cs(a)[ds 2 (a) + ns 2 (a)]cn(x) — ns(a)[cs 2 (a) + ds 2 (a)]cn(x + a) 

— cs(a)ds(a)sn(x + a)dn(x + a) — ds(a)ns(a)sn(x)dn(x) , (174) 

msn(x)cn(x)dn(x + a)sn(x + a) = ds(a)[cs 2 (a) + ns 2 (a)]sn(x) — cs(a)[ds 2 (a) + ns 2 (a)]sn(x + a) 

+ ds(a)ns(a)cn(x + a)dn(x + a) + cs(a)ns(a)cn(x)dn(x) . (175) 

Appendix D: Some Examples of Local Identities of Rank 5 



mdn(x)sn(x)cn(x)dn(x + a)sn(x + a) = — cs(a)ns(a)[cs 2 (a) + ds 2 (a) + ns 2 (a)]cn(x) 



+ 



ns 2 (a)(ds 2 (a) + cs 2 (a)) + cs 2 (a)ds 2 (a) cn(x + a) + cs(a)ds(a)ns(a)sn(x + a)dn(x + a) 



+ds(a)[cs 2 (a) + ns 2 (a)]sn(x)dn(x) + mcs(a)ns(a)cn 3 (x) , 



(176) 



mdn(x)sn(x)cn(x)dn(x + a)cn(x + a) = cs(a)ds(a)[cs 2 (a) + ds 2 (a) + ns 2 (a)]sn(x) 

— ns 2 (a)(ds 2 (a) + cs 2 (a)) + cs 2 (a)ds 2 (a) sn(x + a) + cs(a)ds(a)ns(a)cn(x + a)dn(x + a) 
+ns(a)[cs 2 (a) + ds 2 (a)]cn(x)dn(x) + mcs(a)ds(a)sn 3 (x) , (177) 



m 2 dn(x)sn(x)cn(x)sn(x + a)cn(x + a) = — ds(a)ns(a)[cs 2 (a) + ds 2 (a) + ns 2 (a)]dn(x) 



+ 



ns 2 (a)(ds 2 (a) + cs 2 (a)) + cs 2 (a)ds 2 (a) dn(x + a) + mcs(a)ds(a)ns(a)cn(x + a)sn(x + a) 



+mcs(a)[ds 2 (a) + ns 2 (a)]cn(x)sn(x) + ds(a)ns(a)dn 3 (x) , 



(178) 



mdn(x)sn(x)cn(x)dn 2 (x + a) = cs(a)ds(a)ns(a)[dn 2 (x + a) + 2dn 2 (x) — dn 2 (a) — 2] 

- mcs 2 (a)dn(x)sn(x)cn(x) + [ns 2 (a) (ds 2 (a) + cs 2 (a)) + cs 2 (a)ds 2 (a)] [Z(x + a) - Z(x) - Z(a)][L79) 
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In the following identities, B = — cs 2 (a), B\ = ns 2 (a), B2 = — ds 2 (a). 



dn 2n (x)dn(x + a) = B n dn(x + a) + [ds(a)ns(a)dn(x) - mcs(o)cn(x)sn(x)] ^ B k - l [dxi{x)} 2{n - k ^ , (180) 



fc=i 



dn 2n+1 (x)dn(x + a) = £ n (dn(a) + cs(a)[Z(x + a) - Z(x) - Z(o)]) 

+ [ds(a)ns(a)dn(x) - m cs(a)cn(x)sn(x)] ^ B k ~ l [dn(x)] 2 ( n - fc ) +1 , (181) 



k=i 



m n sn 2n (x)sn(x + a) = B^sn(x + a) - [cs(a)ds(a)sn(x) + ns(a)cn(x)dn(x)] m™"^ 1 [sn(x)] 2(n_fe) , 

(182) 



k=i 



m n+1 sn 2n+1 (x)sn(x + a) = -B?ns(a) [Z(x + a) - Z(x) - Z(a)] 

n 

- [cs(a)ds(a)sn(x) + ns(a)cn(x)dn(x)] ^ m [ sn (x)] 2 ( n - fe ) +1 , (183) 

k=i 

m n cn 2n (x)cn(x + a) = B^c^x + a ) 

n 

+ [cs(a)ns(a)cn(x) - ds(a)sn(x)dn(x)] ^ m^B^-^cn^)] 2 ^-*) , (184) 

k=l 

m n+1 cn 2n+1 (x)cn(x + a) = B 2 n {mcn(a) + ds(o) [Z(x + a) - Z(x) - Z(a)]} 

n 

+ [cs(a)ns(a)cn(x) - ds(a)sn(x)dn(x)] ^ m n ' k B k ' x [cn(x)} 2( - n -^ +1 , (185) 

k=i 

m n cn 2n (x)sn(x)dn(x + a) = -.B^ ns ( a ) cn ( x + a) + cs(a)m n cn 2n+1 (x) 

n 

— ns(a) [cs(a)ns(a)cn(x) — ds(a)sn(x)dn(x)] £ m n - fe B£- 1 [cn(z)] 2(n_fe) , (186) 

fe=i 

m n sn 2n (x)cn(x)dn(x + a) = 5" ds(a)sn(x + a) - cs(a)m n sn 2n+1 (x) 

n 

-ds(o) [cs(a)ds(a)sn(x) + ns(a)cn(x)dn(x)] ]T m n - fc B*- 1 [sn(x)] 2 ( n - fc ) , (187) 

fc=i 

mdn 2n (x)cn(x)sn(x + a) = -B n ds(a)dn(x + a) + ns(a)dn 2n+1 (x) 

n 

-ds(a) [ds(a)ns(a)dn(x) - m cs(a)cn(x)sn(x)] ^ B h_1 [dn(x)] 2 ( n ~ fc ) , (188) 

k=i 
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mdn 2n (x)sn(x + a)cn(x + a) = mB n sn(x + a)cn(x + a) — 2nB n 1 ds(a)cs(a)ns(a)dn(x + a) 

n 

+ ( m + 2ds 2 (a))cs(a) + 2(k - l)ds 2 (a)nc(a)ns(a) 

[ dn ( x )]2(n-fe)+l 



k=l 



— mds(a 

k=i 



(189) 



mdn 2ra+1 (x)sn(x + a)cn(x + a) = — mds(a )ns(a) ^(2fc - l) J B fc - 1 cn(x)sn(a;)[dn(x)] 2 ( n - fe ) +1 



k=i 



+B n - L cs{a) 



(1 - m) - (2n + l)ds(a)ns(a)(dn(a) + cs(a)[Z(x + a) - Z(x) - Z(a)]) 



+ Y,B k - 1 [(m + 2ds 2 (a))cs(a) + 2(fc - l)ds 2 (a)nc(a)ns(a)][dn(a;)] 2 ( n - fc+1 )-S n cs(a)dn 2 (x + 090) 



dn 2n (x)sn(x + a)dn(x + a) = l? n sn(x + a)dn(x + a) - 2ni? n 1 ds(a)cs(a)ns(a)cn(x + a) 

n 

- Z! B * _1 [ cs ( a ) ns («) + 2(fc - l)ds 2 (a)nc(a)] sn(x)[dn(x)] 2 ( n - fe ) +1 
k=l 

n 

+ds(a) ^^ 1 [cs 2 (a) + (2A; - l)ns 2 (a)]cn(x)[dn(x)] 2 ( n - fe ) , 



(191) 



k=i 



dn 2n+1 (x)sn(x + a)dn(x + a) = -B n cs(a)cn(x + a)dn(x + a) + £ n_1 cs(a)ns(a)[cs 2 (a) + 2nds 2 (a)]sn(x) 

n 

- ^ J B fc - 1 [cs(a)ns(a) + 2(fc - l)ds 2 (a)nc(a)]sn(x)[dn(x)] 2 ( n " fc+1 ) + (2n + l)5 n ds(a) ns(a)sn(x + a) 



fc=i 



+ds(a) ]T 5 fc - 1 [cs 2 (a) + (2A; - l)ns 2 (a)]cn(x)[dn(x)] 2(n - fe)+1 , 



(192) 



k=i 



dn 2n (x)cn(x + a)dn(x + a) = B n cn(x + a)dn(x + a) + 2n£? n 1 ds(a)cs(o)ns(a)sn(x + a) 

n 

- B^ 1 [cs(a)ds(a) + 2(k - l)ds(a)nc(a)ns(o)] cn(x)[dn(x)] 2 ( n - fc ) +1 



k=i 



-ns(a) J] ^ fc_1 [cs 2 (a) + (2k - l)ds 2 (a)]sn(x)[dn(x)] 2 ( n - fc ) 
k=i 



(193) 



dn 2n+1 (x)cn(x + a)dn(x + a) = 5 n cs(a)sn(x + a)dn(x + a) + £ n_1 cs(a)ds(a)[cs 2 (a) + 2nns 2 (a)]cn(x) 

n 

- 5] J B fe - 1 [cs(a)ds(a) + 2(k - l)ds(a)nc(a)ns(a)]cn(x)[dn(x)] 2(n - fe+1) + (2n + l)£ n ds(a) ns(a)cn(x + 



fc=i 



-ns(a) ^ J B fc - 1 [cs 2 (a) + (2ife - l)ds 2 (a)]sn(x)[dn(x)] 2(n - fe)+1 
k=i 



(194) 
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dn 2n (a;)dn 2 (x + a) = 2n£"- 1 ds(a)ns(a)[dn(a) + cs(a)(Z(x + a) - Z(x) - Z(a))] 

n-l 

-(1 - m)^- 1 + B k -\B 2 - (1 - m) + 2Ms 2 (a)ns 2 (a)][dn(x)] 2 ( n - fc ) 



fe=i 

n-l 



+ J B n dn 2 (x + a) + £dn 2n - 2mcs(a)ds(a)ns(a)sn(a;)cn(x) ^ feS^^dn^)] 2 ^-^" 1 , (195) 



k=i 



dn 2n+1 (x)dn 2 (x + a) = Bdn 2n+1 (a;) + B n mcs(a)cn(x + a)sn(x + a) 

n 

+ ^ S^fS 2 - (1 - m) + 2Ms 2 (a)ns 2 (a)][dn(x)] 2(n - fc)+1 



k=i 

n 

— 2mcs(a)ds(a)ns(a)sn(x)cn(x) kB^idnix)] 2 ^-^ + (2n + l)£ n ds(a)ns(a)dn(x + a) , (196) 

fc=i 

Appendix F: Examples of Identities with Weighted Terms and Their 
Linear Combinations 

In this appendix a = ^f, a' = b = ^f. 
Fl: MI-I Identities 



p v 
m S c 'A s i+r ~ s j-r] = 2[ns(a) - ds(a)] 53 4/ > ( 197 ) 
j"=i j"=i 
p p 

d|[dj+ r — rfj-r] = — 2mcs(a) 53 CjSj , (198) 

p p 

53 Cj [c j+r d j+r - Cj-rdj-r] = 2ds(a) 53 c jSj , (199) 

3=1 3=1 
V P 

Sj [sj +r dj +r — Sj- r dj- r ] = — 2ns(a) 53 CjSj , (200) 

3=1 3=1 



p 

J] ^ [d j+r c j+s - dj- r Cj- s ] = , (201) 

3=1 
V 

53 Sj[d j+r s j+s - dj-rSj- a ] = , (202) 

3=1 
P 

53 dj [d j+r d j+s - dj- r dj- s } = , (203) 

3=1 
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dj [ c j+rCj+s Cj-rCj-s] — , 

3=1 
P 

^ ' dj[sj+ r Sj+ s Sj— r Sj— s ] = , 

3=1 

V V 

m ^ dj[cj +r Sj +r — Cj^ r Sj^ r ] = 2cs(a) [ns 2 (a) + ds 2 (a) — 2ns(a)ds(a)] ^ dj , 



3=1 



3=1 



m ^ Cjdj[dj +r Sj +r — dj- r Sj- r ] = — 2[ns(o)(cs 2 (a) + ds 2 (a)) — ds(a)(cs 2 (a) + ns 2 (a))] ^ dj , 

V V 

d| [d| +r — d|_ r ] = — 2mcs(a)[cs 2 (a) + 2ds(a)ns(a)] ^ CjSj , 
i=i i=i 

p v 
m ^ CjSjdj[cj +r Sj +r — Cj^ r Sj^ r ] = 2cs(a)[(ds 2 (a) + ns 2 (a)) + ds(a)ns(a)] CjSj , 

3=1 3=1 



(204) 
(205) 
(206) 
(207) 
(208) 
(209) 



mJ2cjSjdj[dj+ r ~ dj-r] = 2cs(a) ^ d) 

3=1 3=1 

V 

-2cs(a)[ds 2 (a) + 1 - ds(a)ns(a)] ^ dj , 

3=1 



(210) 



V V 

m ^ c j s j4/[^j+r ~~ ^|-r] = 2cs(a)[ds(a)ns(a) — cs 2 (a)] ^ d| + 2cs(a) 

3=1 ' 3=1 



2cs 2 (a)ds 2 (a 



+2ns 2 (a)(ds 2 (a) + cs 2 (a)) - cs 4 (a) - ds(a)ns(a)(ds 2 (a) + cs^(a) + ns 2 (a)) 



(211) 



m Cjdj[dj +r Sj +s — (fj_ r Sj_ s ] = 2[ns(a — a')cs(a)(ds(a) — ns(a)) — cs(a — a')cs(a')(ds(a') — ns(a'))] Y^ dj , 

3=1 3=1 

(212) 

P P 

m Y~^ Sjdj[dj +r Cj +s — dj- r Cj- s ] = — 2[ds(a — a')cs(a)(ds(a) — ns(a)) — cs(a — a')cs(a')(ds(a / ) — ns(a'))] ^ dj , 

(213) 

p p 
m ^2 CjSj[dj+ r dj+ s — dj^ r dj- s ] = 2cs(a — a')[ns(a)(ds(a) — ns(a)) — ns(a')(ds(a') — ns(a'))] ^ dj , (214) 
j=i j"=i 
p p 
m 2 ^ CjSj[sj +r Sj +s — Sj_ r Sj_ s ] = — 2ns(a — a')[ns(a)(ds(a) — ns(a)) — ns(a')(ds(a') — ns(a'))] ^ dj , (215) 

3=1 3=1 

p p 
m 2 Y~^ CjSj[cj +r Cj +s — Cj- r Cj- s ] = 2ds(a — a')[ns(a)(ds(a) — ns(a)) — ns(a')(ds(a') — ns(a'))] Y~] dj , (216) 
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P V V 

dj[dj+ r — dj^ r ] = — 2mcs(o) ^ c j s j^j + 2mcs 3 (a) ^ CjSj , (217) 
j=i j=i j=i 

VP V 

m 2 c)[s j+r - Sj -r] = 2ns(a) ^ d) - 2[ns(a)(2ds 2 (a) - cs 2 (a)) - ds 3 (a)] ^ , (218) 
j=i j=i j=i 



v 



m 2 J2 s^Cj+r - Cj- r ] = -2ds(a) ^ 4 + 2 [ds(a)(ns 2 (a) + 1) - ns 3 (a)] J2 d j > ( 219 ) 



F2: MI-II Identities 



p p 
m ^2 Sjdj[cj +r — Cj- r ] = 2ds(a) ^ d 2 — 2pns(a)[dn(a) — cs(a)Z(a)] , (220) 
i=i 3=1 



v v 



m ^ Cjdj[sj +r — Sj- r ] = 2ns(a) ^ d 2 - — 2pds(a)[dn(a) — cs(a)Z(a)] , (221) 

3=1 3=1 
V V 

m ^2 SjCj[dj+ r — dj- r ] = 2cs(a) ^ dj — 2p[cs(a) — ds(a)ns(a)Z(a)] , (222) 

3=1 3=1 

V 

m ^2 dj [sj +r Cj +s — Sj_ r Cj_ s ] = — 2pds(a — a')[dn(a) — cs(a)Z(a)]+2pns(a — a')[dn(a') — cs(a')Z(a')] , (223) 

3=1 
V 

m Cj[(ij +r Sj +s — dj- r Sj- s ] = 2pdc(a)ns(a— a')[dn(a)— cs(a)Z(a)]— 2pdc(a')cs(a— a')[dn(a') — cs(a')Z(a')] , 
i=i 

(224) 

p 

m ^ Sjfdj+rCj+s— dj- r Cj- s ] = 2pnc(a)ds(o— a')[dn(a)— cs(a)Z(a)]— 2pnc(a')cs(a— a')[dn(a') — cs(a')Z(a')] , 
i=i 

(225) 

p p 

d 3 [dj +r — dj^ r ] = — 2mcs(a) ^ dj s j c j j (226) 

p p 

111 Z 4 \- S 3+r ~ S 3~r] = -2ns(a) J2 d 3 S 3 C 3 > ( 227 ) 
P P 

m ^ c 3 [cj +r — Cj_ r ] = — 2ds(a) ^ djSjCj , (228) 
i=i i=i 



m ^ d 2 [dj-|_ r s :) ' +r c :) ' +r — dj- r Sj_ r Cj- r ] = — 6cs(a)ds(a)ns(a) ^ d 2 



p p 

i 2 

i=i 3=i 
+2pcs(a)ds(a)ns(a)[dn 2 (a) + 2] - 2p[ns 2 (a)(cs 2 (a) + ds 2 (a)) + cs 2 (a)ds 2 (a)]Z(a) , (229) 
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p 

m SjCj[dj +r — dj_ r \ = — zcs(a)[ds (a) + ns (a) + 


p 

cs (a) J 2^ dj 

3=1 




-2pcs(a)[l — m — 3ds 2 (a) + 3cs(a)ds(a)ns(a)Z(a)] , 




(230) 


p 

m "i c ii s j+r ~~ s j-r\ = 2ns(a)|ds (a) + ns (a) + < 


V 

3S (a)J 2^ a,- 

3=1 




+2pns(a)[l — m — 3ds 2 (a) + 3cs(a)ds(a)ns(a)Z(a)] , 




(231) 


p 

m 2 ^ djSj[cf +r - c?_ r ] = -2ds(a)[ds 2 (a) + ns 2 (a) + 
i=i 


V 

cs 2 (a)]£d 2 




— 2ods(a)[l — m — 3ds 2 (a) + 3cs(a)ds(a)ns(a)Z(a)l , 

■L \/L \ / \/\/\/\/J? 




(232) 


?3: MI-III Identities 






P V 

dj[c j+r - Cj _ r ] = 2[cs(6) - ds(6)] ^ Sj , 




(233) 


p p 
^ ms 2 [s j+r - sj-_ r ] = -2ns(6) ^ Cj-dj , 




(234) 


p 

^ ' Cj [Sj+ r Cj+ s — Sj_ r Cj_ s ] = , 

3=1 




(235) 


p 

^ ] (ij [<ij_|_ r Sj_|_ s — dj_ r Sj_ s ] = , 
3=1 




(236) 


p 

s j[dj+rdj+s — dj- r dj- s ] = , 




(237) 


p 




(238) 


p 

s j [ s j+rSj+s — Sj_ r Sj_ s ] = , 

3=1 




(239) 


P P 

(^[cj+rdj+r - Cj- r dj- r } = -2ns(6) [cs 2 (6) + ds 2 (6) - 2cs(6)ds(6)] Sj , 

j=l 3=1 


(240) 


V 

m ^ Sjdj [cj +r Sj +r — Cj- r Sj- r ] = — 2[ds(6)(cs 2 (6) + ns 2 (6)) — cs(6)(ds 2 

3=1 


(6)+ns 2 (&))]]>> , 

3=1 


(241) 


V V 

m 2 J2 s M+r - s j-r] = 2ns(6)[ns 2 (6) + 2ds(6)cs(6)] ^ Cjdj , 


(242) 



3=1 3=1 



P P 

m ^2i c j s jdj[ c j+rdj+r — Cj- r dj- r ] = 2ns(6) [(ds 2 (b) + cs 2 (6)) + ds(6)cs(6)] V] Cjdj , (243) 

3=1 3=1 
P P 

mJ2 c j s j d j[ s j+r - Sj-r] = -2mns(6) »? 

3=1 3=1 

P 

-2ns(6)[(-ds 2 (&) + 1)) + ds(6)ns(6)] ^ Sj , (244) 

3=1 



m2 c 3 s J d A s3 j+r ~ s i-r] = 2m ns(6) [ds(6)cs(fe) - ns 2 (6)] ^ s 3 - 2ns(6) 
j=i i=i 

+ns 2 (6)(2ds 2 (6) + cs 2 (6)) - 2ds(6)cs(6)(ds 2 (6) + cs 2 (6) + ns 2 (6)) 



3cs 2 (6)ds 2 (6) 



J2 S 3 > ( 245 ) 



m ^ CjSj^j+rSj+s — (ij_ r Sj_ s ] = 2[ns(6 — 6')ns(6)(ds(6) — cs(6)) — cs(6 — 6')ns(6')(ds(6 / ) — cs(6'))] ^ Sj , 

3=1 3=1 

(246) 

p p 
m ^2 Sjdj[sj+ r Cj+ s — Sj-rCj-s] = — 2[ds(6 — 6')ns(6)(ds(6) — cs(6)) — ns(6 — 6')ns(6')(ds(6 / ) — cs(6'))] ^ Sj , 

3=1 3=1 

(247) 

p p 
J^Cjdjldj+rdj+s - dj-rdjs] = -2cs(6 - b')[cs(b)(ds(b) - cs(6)) - cs(&')(ds(6') - cs(6'))] Yl s i > ( 248 ) 

3=1 3=1 
P P 

m Cjdj[s j+r s j+s - Sj-rSjs] = -2ns(6 - b')[ds(b)(ds(b) - cs(b)) - ds(6')(ds(6') - cs(6'))] s i > ( 249 ) 

3=1 3=1 

p p 
m J2 Cjdj[c j+r c j+s - Cj-rCjs] = 2ds(6 - b')[ds(b)(6s(b) - cs{b)) - ds(6')(ds(6') - cs(b'))] ^ Sj , (250) 

3=1 3=1 

p p p 

m 2 J2 sfaj+r - Sj-r] = 2ns(6) c 3 d ) - 2ns(6)[ns 2 (6) + 1] c j d j , (251) 

3=1 3=1 3=1 

P P P 

d)[c j+r - Cj- r ] = 2mds(6) ^ s| — 2[cs 3 (6) - ds(6)(ns 2 (6) - 2] ^ Sj , (252) 

3=1 3=1 3=1 

P P P 

m c )[dj+r ~ dj-r] = 2mcs(6) ^ s 3 - 2[ds 3 (6) - cs(6)(2ds 2 (6) - 1] ^ Sj , (253) 

3=1 3=1 3=1 

F4: MI-IV Identities 



p p 

dj [s j+r ~ Sj-r] = 2[ns(6) - cs(6)] ^ cj , (254) 

3=1 3=1 
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P P 

m c 2 [c ? +r — c j-r] = — 2ds(6) ^2 s jdj , (255) 
j"=i 3=1 
p 

^2 s j[ s j+r c j+s ~ Sj-rCj-s] = , (256) 
3=1 
P 

d j [dj+rCj+s ~ dj-rCj-s] = , (257) 

3=1 

p 

^2 Cj[dj+ r dj +s - dj- r dj- s ] = , (258) 

3=1 

p 

c j [ c j+r c j+s ~ c j-r c j— s] = , (259) 

3=1 

p 

^2 c j[ s j+r s j+s — s j-r s j-s] = , (260) 
3=1 

P P 

J2 djisj+rdj+r - sj^rdj-r] = 2ds(6) [cs 2 (6) + ns 2 (6) - 2cs(6)ns(6)] J2 c j > ( 261 ) 

3=1 3=1 

p p 
m ^2 Cjdj[cj +r Sj +r — Cj- r Sj- r } = 2[cs(6)(ds 2 (6) + ns 2 (6)) — ns(6)(ds 2 (6) + cs 2 (6))] ^ Cj , (262) 

3=1 3=1 
P P 

m 2 J2 c 3 j[£j+r - c 2 j~ r ] = -2ds(6)[ds 2 (6) + 2ns(6)cs(6)] J2 Sj dj , (263) 

3=1 3=1 

p p 
m^2 CjSjdj[sj +r dj +r — Sj^ r dj- r ] = 2ds(6)[(ns 2 (6) + cs 2 (6)) + ns(6)cs(6)] V] Sjdj , (264) 

3=1 3=1 



p 

3=1 3=1 



m J2 c jSjdj[cj+r - Cj_ r ] = 2mds(6) ^2 c j 



p 



-2ds(6)[(cs 2 (&) + m)) - cs(6)ns(6)] J2 c i > ( 265 ) 

3=1 



p p 

2/u\ i „„2/ 



m2 H 'v-syJ, (••'. ,. - c|_ r ] = 2mds(6)[ns(6)cs(6) - ds 2 (6)] ^ c i ~ 2ds ( & ) 
j'=i j"=i 

p 



2cs(b)ns(b)(ds 2 (b) + cs 2 (b) 



+ns 2 (b)) - ((2ds 2 (6) + 3cs 2 (6))ns 2 (6) + ds 2 (6)(2cs 2 (6) - ds 2 (b)) 



J2c 3 , (266) 
j'=i 



CjSj[dj +r Cj +s — dj- r Cj- s ] = 2[ds(6 — 6')ds(6)(ns(6) — cs(b)) — cs(b — b')ds(b')(ns(b') — cs(6'))] ^ Cj , 

3=1 

(267) 
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P P 

m Cjdjlsj+rCj+s - Sj-rCjs] = -2[ds(6 - 6')ds(6)(ns(6) - cs(6)) - ns(6 - b')ds(b')(ns(b') - cs(b'))] ^ cj , 

i=i 3=1 

(268) 

p p 
J2 djSj[d j+r d j+s - dj-rdjs] = 2cs(6 - b')[cs(b)(ns(b) - cs(6)) - cs(6')(ns(&') - cs(6'))] E C J > ( 269 ) 
i=i i=i 
p p 
m ^2 djSj[sj +r Sj +s — Sj- r Sj-. s ] = 2ns(6 — 6')[ns(6)(ns(6) — cs(b)) — ns(6')(ns(6') — cs(6'))] ^ Cj , (270) 

3=1 3=1 
P P 

m Y djSj[c j+r c j+s - Cj-rCjs] = -2ds(6 - 6')[ns(6)(ns(6) - cs(6)) - ns(b')(us(b') - cs(b'))] E cj , (271) 

3=1 3=1 
P P P 

m 2 Cj[cj+r - cj-r] = -2ds(6) E s 3 d ) + 2ds(6) [2ds 2 (6) - cs 2 (fe)] ^ , (272) 

3=1 3=1 3=1 

P P P 

E d)[s j+r - Sj - r ] = 2mns(6) ^ cf + 2[cs 3 (6) + ns(6)(ds 2 (6) - 2cs 2 (6)] ^ Cj , (273) 

3=1 3=1 3=1 

P P P 

m s )[dj+r ~ dj-r] = -2mcs(6) ^c 3 - 2[ns 3 (6) - cs(6)(ns 2 (6) + m] ^ Cj , (274) 
j=i j=i j=i 

F5: MI-I Identities With Alternating Signs 

Note in this and the next subsection, p is even and r is necessarily an odd integr. 

p p 
m E(- 1 )'" ls i[ c i+- " c 3-r\ = 2 ( ds («) + ns («)) > ( 275 ) 

3=1 3=1 

j2(-iy~ l d)[d 3+r - dj- r ] = -2mcs(a) £(-1)'- , (276) 

3=1 3=1 

j^i-iy^cjdjici+r - cj- r ] = -2ds(a) £(-1)'- , (277) 
E(-l) J ' _1 «idi[si+r " sj-r] = 2ns(a) X^-l)'" 1 ^ , (278) 

3=1 3=1 
P P 

Y(- 1 ) i ~ 1 d 2 j [cj+ r s j+ r - Cj-rSj-r] = 2cs(a) [ds 2 (a) + ns 2 (a) + 2ds(a)ns(a)] 5^(-l) J ' -1 dj , (279) 

3=1 3=1 

p p 
E( — Sjdj[cj +r dj +r — Cj- r dj- r ] = 2[ns(a)(cs 2 (a) + ds 2 (a) — ds(a)(cs 2 (a) + ns 2 (a)] E( — ^y~ lc j s j > 

3=1 3=1 

(280) 

F6: MI-II Identities with Alternating Signs 
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p p 
m ]T(-ir V^+r - d,_ r ] = 2cs(a) ^(-lr 1 ^ , 

3=1 3=1 




(281) 




P V 

m ^-ly-^djSjicj+r - Cj . r ] = 2ds(a)^(-ir 1 4 , 

J'=l j=l 




(282) 




p p 
m £(-irV> j+r - aj _ r ] = 2ns(a) ^-l)'"^ , 

3=1 3=1 




(283) 




P P 

m ]T( W ■r J s J d). r - 4_ r ] = 2cs(a)(ds 2 (a) + 1) £(-1)'- ^ , 

3=1 3=1 




(284) 




P P 

m ^(_i),-i d . s .[ c 3 +r _ c 3_ r] =2 ds(a)(cs 2 (a)+m)^(-ir 1 ^ , 
j"=i j"=i 




(285) 




m £(-1)'- V>?+r - *?-r] = M«)(l " dg2 («)) E(-1) J_1 4 > 
J"=l J"=l 




(28b) 


V P V 

5^(-l) J ' -1 <§[d j+r - dj-r] = -2mcs(a) ^(-ly^djSjCj + 4cs 3 (a) 
j=i j=i j=i 


i-l 7 


(287) 


p 

2 E(- 


P P 

-l) J ' -1 s|[sj +r - Sj- r ] = -2mns(a) ^(-l) J ' _1 djSjCj + 4ns 3 (a) 


-17. 

3 ' 


(288) 


2 E(- 


-ly-^cj+r-cj-r] = -2mds(a)£(-iy- 1 ^ SjCj +4ds 3 (a)X:(-l) r 

j"=i j"=i 

^(-ly- 1 ^ 3 ^ - 4_ r ] = 4mcs 3 (a) f^-l)'- 1 ^- 
i=i j=i 


3 ' 


(289) 




p 

-4cs(a)[ns 2 (a)(cs 2 (a) + 3ds 2 (a)) + cs 2 (a)(cs 2 (a) + ds 2 (a))] ]T(-l) j ~ 

i=i 




(290) 




m 3 E(-iy- 1 4[4 +r - S 3 _ r ] = -4mn S 3 (a)£(-iy- 1 Q !, SjCj 

3=1 3=1 








P 

+4ns(a)[ns 2 (a)(cs 2 (a) + ds 2 (a) + ns 2 (a)) + 3cs 2 (a)ds 2 (a)] ^(-l) j_] 

j'=i 




(291) 




m 3 El-iy-'c 3 ^ - c 3 _ r ] = 4mds 3 (a) £( 1H %s jCj 

3=1 3=1 








P 

-4ds(a)[ns 2 (a)(3cs 2 (a) + ds 2 (a)) + cs 2 (a)(cs 2 (a) + ds 2 (a))] 

3=1 




(292) 
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m 2 ^(— iy l djSjCj[dj +r Sj +r Cjj rr — dj- r Sj- r Cj_ r ] = — 4mds(a)cs(a)ns(a) y^(— 1) J l djSjCj 

3=1 3=1 

V 

+8ds(a)cs(a)ns(a)[ns 2 (a) + cs 2 (a) + ds 2 (a)] ^(-ly^Zj , (293) 

3=1 

m 2 ^(_iy-i d . s 2 c . [s . +r _ s ._ r] = _ 2ns(a) + 2ns(a)[ds 2 (a) + 1] jj-l^^ , (294) 

j'=i i=i 3=1 

Appendix G: Some Definite Integrals 

Here K, E correspond to the complete elliptic integrals of the first and second kind respectively. 

,-2K 

/ dn 3 (a;)dn(a; + a)dx = 2ds(a)ns(a) J B-2 J ft:cs 2 (a)[dn(a)-cs(a)Z(a)], (295) 
Jo 

,-2K 

/ m 2 sn 3 (x)sn(x + a)dx = 2cs(a)ds(a)E - 2K[cs(a)ds(a) - ns 3 (a)Z(a)], (296) 
Jo 

,-2K 

/ m 2 cn 3 (x)cn(x + a) dx = 2cs(a)ns(a)£ + 2K[m 2 cn(a) - cs(a)ns(a) + ds 3 (a)Z(a)], (297) 
Jo 

,-2K 

I mdn(x)sn(x)dn(x + a)sn(x + a) dx = 4cs(a)ns(a)_E — 2ATns(a)[cs(a)(l + dn 2 (a)) — (1 + cn 2 (o))ds(a)ns(a)Z(o)], 
Jo 

(298) 

/ mdn(x)cn(x)dn(x + a)cn(x + a) dx = -4cs(a)ds(a)E + 2K[2cs(a)ds(a) - (cs 2 (a) + ds 2 (a))ns(a)Z(a)], (299) 
Jo 

r-2K 

/ m 2 sn(x)cn(x)sn(x + a)cn(x + a) dx = 4ds(a)ns(a)-E + 2ATns(a)(l + dn 2 (a))[cs(a)ns(a)Z(a) — ds(a)], (300) 
Jo 

1 /" 2K 

— — / dn(x)dn(x + a)dn(x + a )dn(x + a") dx = dn(a)dn(a')dn(a") + cs(a)cs(a' — a)cs(a" — a)Z(a) 
2 a Jo 



1 /" 

— / m 2 sn(x)sn(x + a)sn(x + a')sn(x + a") = ns(a)ns(a' — a)ns(a" — a)Z(a) 
Jo 

- ns(a')ns(a' - a)ns(a" - a')Z(a') + ns(a")ns(a" - a)ns(a" - a')Z{a"), (302) 



1 

2^ 



2K 

m 2 cn(x)cn(x + a)cn(x + a')cn(x + a") g?x = m 2 cn(a)cn(a')cn(a") + ds(a)ds(a' — a)ds(a" — a)Z(a) 



- ds(a')ds(a' - a)ds(a" - a')Z(a') + ds(a")ds(a" - a)ds(a" - o')Z(o"), (303) 

1 /- 2K 

— — : / m 2 cn(x)sn(x + a)cn(x + a')sn(:r + a") = m 2 sn(a)cn(a')sn(a") — ds(a)ds(a' — a)ns(a" — a)Z(a) 
Jo 

+ ds(a')ns(a' - a)ns(a" - o')Z(o') - ds(a")ns(a" - a)ds(a" - a')Z(a"), (304) 
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2 A' 



/ msn(x)dn(i + a)sn(x + a')dn(x + a") dx = — ns(a)ns(a' — a)cs(a" — a)Z(a) 

2K Jo 

+ ns(a )cs(a — a)cs(a" — a')Z(a') — ns(a")cs(a" — a)ns(a — a )Z(a ), (305) 



1 /- 2if 

— — / mcn(i)dn(s + a)cn(x + a )dn(x + a") dx = nidn(a)cn(a')dn(a") + ds(a)ds(a' — a)cs(a" — a)Z(a) 



ds(a')cs(a' - a)cs(a" - a')Z(a') + ds(a")cs(a" - a)ds(a" - a')Z(a"). (306) 



Appendix H: Some Indefinite Integrals 

We give below recursion relations expressing certain arbitrary order integrals in terms of lower order 
integrals, well known integrals of sn n (x), dn n (x), cn n (x) [3], and incomplete elliptic integrals of the first, 
second and third kind, which essentially occur due to the integral (|67|). It should be noted that in view 
of the identities (|88[) and (|89|). the integrals for msn(x + a)sn(x) and mcn(x + a)cn(x) are related to the 
integral (|67|). In this appendix, n > 1. 

7» = ns 2 (a)J n _i - m " 1dS ^ sn 2 "-^!) - cs(a)ds(a)m"- 1 f sn 2n - 1 (x)dx, (307) 
(2n-l) 7 

where, Ik = J m' c sn 2fe (x)sn(a; + a) (ice. 

J w = -ds 2 (a)/»_i + m " lds( ° ) cn 2 "- 1 (x)+cs(a)ns(a)m"- 1 / cn 2 "- 1 ^) dx, (308) 
(2n-l) 7 

where, 7^ = J" m fe cn 2fc (x)cii(x + a) da;. 

7„ = ns 2 (a)7 n _ a - mn ~' ds ( a )° s ( a ) sn 2 "- 1 ( 3; ) + cs(a)m n / sn 2 "- 1 ^ 2 ^) dx, (309) 



(2ti - 1) 

where, Ik = J m sn (x)cn(a;)dii(a; + a) dx, 7q = ds(o) f sn(x + a) dx — cs(a) f sn(x) dx 



7„ = -ds 2 (a)/„_ 1 - m " lds ( a )^ s (a) cn 2n-i (x) _ cs(a)m n / cn 2 "- 1 (x)sn 2 (x)dx, (310) 



(2n - 1) 

where, Ik = J na cn 2fc (x)sn(x)dn(x + a) dx, Iq = — ns(a) J" cn(x + a) dx + cs(a) / cn(x) dx 



I n = -cs 2 (a)J n _i - CS } a)dS ^ dn 2n - 1 (x)+mns(a) j dir" 'i.rinr(..-)d.r. (:5II) 
(2n — 1) 
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where, h = J mdn 2fe (x)cn(x)sn(x + a) dx , Jo = — ds(a) J dn(x + a) dx + ns(a) / dn(x) dx. 

I n = BI^ + d ^ a ) ns ( a ) dn 2n-i (a; ) + 2 ds(a)ns(a) V [dii(a?)] 2 t" fc ) _ 2cg(a)ds(a)ns(a) f dn , x + a) dx 
(2n-l) J 2(n-fc)-l 7 

/n— 1 « 
dn 2 "" 1 ^) dx - 2cs(a)ds 2 (a)ns 2 (a) ^ B^ 1 / [dn(aO] 2(n_fc)_1 dx, (312) 
/, — i 



where, Ik = / mdn (x)cn(x + a)sn(x + a) dx, £> = — cs 2 (a) , J = — dn(x + a). 

/„ = -cs 2 (a)J„_i + £^dn 2 "(x) + ds(a)ns(a) / dn 2 "(x) dx , (313) 
where, h = J mdn 2fe+1 (x)dn(x + a) dx , lo = J dn(x + a)dn(x) dx . 

J„ = ns 2 (a)/„_i - m " ( 2° ) S(fl) sn2 "( x ) - m"- 1 cs(a)ds(a) J sn 2 "(x) dx , (314) 
where, Ik = J ni fe+1 sn 2fe+1 (x)sn(x + a) dx , I = J msn(x + a)sn(x) dx = +dc(a)x — nc(a) / dn(x + a)dn(x) dx. 

/„ = -ds 2 (a)/„_i + m " lds(fl) cn 2 "(x) + m"- 1 cs(a)ns(a) / cn 2 "(x) dx , (315) 
where, Ik = Jm fe+1 cn 2fe+1 (x)cn(x + a) dx , 7 = Jmcn(x + a)cn(x) dx = — (1 — m)nc(a)x + ds(a) Jdn(x + a)dn(x) dx. 

/„ = B/„_i + ds(a)n " (a) dn 2 "(x) + ds(a)ns(a) V i? fc f^)] 2 ^ _ 2 c S ( a )ds 2 (a)ns 2 ( a ) V S^ 1 / [dn(x)] 2 ("- fc > dx 
(2n) ^ (n-fc) ^ y 

+cs(a)[m + 2ds 2 (a)] J dn 2 "(x) dx - 2cs(a)ds(a)ns(a)B™- 1 J dn(x + a)dn(x) dx , (316) 
where, h =J nidn 2fe+1 (x)sn(x + a)cn(x + a) dx , £> = — cs 2 (a) , 

Iq = /mdn(x)sn(x + a)cn(x + a)dx = —(1 — m )cs(a)x + cs(a)_E(am x, fc) + ds(a)nc(a) / dn(x + a)dn(x)dx. 
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